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INTRODUCTION 

This unit differs from the others in the course in that it seeks to show how 
complex analysis may be used in the mathematical modelling of a physical 
process, namely, the flow of a fluid such as water or air. 

If you have had experience of mathematical modelhng, then the approach 
adopted here will be familiar to you, in that an initial simpliflcation of a 
real-world situation is required before we axe able to express what is going on 
mathematically and make progress by means of posing and solving problems. 
Later we interpret the solutions in physical terms, and conclude by analysing 
how good the model is and where it could be improved. 

If you have not done mathematical modelling before, you need not feel 
apprehensive about it. This unit may be treated primarily as a mathematics 
text, although it will help your appreciation of the subject matter if you axe 
able to think in terms of physical quantities such as velocity and pressure. 

When mathematical modelling is carried out ‘for real’, it is not always clear at 
the outset which body of mathematical knowledge to apply to the real-world 
problem being faced. The topic of fluid flows, however, has been chosen 
precisely because, when various modelhng assumptions axe made, it can be 
viewed as an apphcation of complex analysis. 

Although the assumptions required to achieve our mathematical end axe quite 
restrictive, the resulting model is stiU significant, both historically and in terms 
of the physical insights which it provides. More ‘realistic’ models of fluid flows 
have been developed, but their additional realism is paid for by an increase in 
mathematical complexity. The model considered here has the virtue of being 
non-trivial but of sufficient simplicity for its predictions to be ihuminating. For 
example, the model predicts correctly the swerve on a spinning ball projected 
through the air, and the presence of the upward force which keeps an aeroplane 
in flight. 

This same simple model of fluid flows can be developed mathematically 
without the benefit of complex analysis, but the natural relationships within it 
cannot then be exhibited in quite such a transparent manner. Complex 
analysis is a natural language with which to describe this particular model. 

In Section 1 we set up the mathematical model for fluid flows, using a 
complex-valued function to describe the steady fluid velocity within a 
cross-sectional plane. It turns out that, under suitable modelling assumptions, 
the conjugate of the velocity function is analytic, this opens up a range of 
possibilities for the apphcation of previously estabhshed mathematical results. 

In particular, the conjugate velocity function has a primitive (on a suitable 
region) known as a complex potential function, which we introduce in Section 2. 
The streamlines, or paths of points moving with the fluid, can be described 
simply in terms of the complex potential. In the audio tape in Subsection 2.2 
we look at several examples of simple fluid flows, and derive their complex 
potentials and the corresponding streamhne patterns. 

In Section 3 we turn to the investigation of flow patterns in a stream passing 
axound a sohd object or obstacle, covering in detail the case of flow past a 
circular cylinder. We then show how the theory of conformal mappings may be 
used to link the results for a circular cylinder to those for obstacles of other 
shapes. 

In certain circumstances, the pattern of flow velocity past an obstacle gives rise 
to a force on the obstacle, and Section 4 demonstrates how this force may be 
calculated. The study of flow past aerofoils, which was introduced as an 
important special case in Section 3, is continued in Section 4. Finally, we talce 
stock of how well the predictions of the model compare with reality. 
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study guide 

Section 3 forms the most important as well as the heaviest part of this unit, 
and you will probably need to spend almost half of your study time on it. You 
should fin d that Section 2 (which contains the audio tape) is relatively light, 
although essential in terms of building up your ability to visualize how fluids in 
motion behave according to the model. Subsections 4.3 and 4.4 axe intended 
for reading only. 


1 SEHING UP THE MODEL 

After working through this section, you should be able to: 

(a) appreciate the modelling assumptions which are made in order to represent 
fluid flow velocity by a continuous complex function; 

(b) apply the formula for the components of velocity at any point, in the 
direction specified by e'®; 

(c) explain what is meant by stagnation point, streamline, unit-speed 
parametrization, circulation, flux, locally circulation-free, locally flux-free, 
model flow, source, sink and vortex’, 

(d) understand why the conjugate of a model flow velocity function is an 
analytic function, and conversely; 

(e) establish whether a velocity function is locally circulation-free, locally 
flux-free, or both. 


1.1 A complex-valued velocity function 

The aim of this unit is to develop a mathematical model which can be applied 
to the flow (or motion) of a fluid. Webster’s Dictionary defines a ‘fluid’ as ‘a 
substance that alters its shape in response to any force however small, that 
tends to flow or to conform to the outline of its container, and that includes 
gases and liquids ... ’. Both water (a liquid) and air (a gas) are fluids, and your 
everyday experience of how these substainces behave in motion should give you 
some feeUng for what we intend to model. 

One way of visualizing a fluid flow is shown in Figure 1.1 below. The diagram 
represents a cross-section, at a particular instant in time, of the flow of some 
fluid (water, say). 



Each arrow in this diagram represents the instantaneous velocity of the fluid at 
the point from which the arrow is directed. The magnitude of this velocity 
(that is, the fluid speed) is represented by the length of the arrow, while the 
flow direction is given by the direction of the arrow. 


Velocity is a vector queintity; 
that is, it has both magnitude 
and direction. 



Photographs of flows rather similar to Figure 1.1 (but without the arrowheads) 
may be obtained by the insertion of small neutrally buoyant beads into the 
flow. Exposure of the film for a short interval of time results in the movements 
of the beads showing up as streaks on the photograph. Since the beads are 
carried along by the fluid, the streak lengths and directions indicate fluid 
velocity, as represented above by the arrows. 

Suppose that we regard the flow cross-section shown in Figure 1.1 as lying in a 
region of the complex plane. Then each point in the region can be described by 
a complex number z = x->riy. Since the flow velocity at the point z, like the 
arrow which represents it, is characterized completely by its ma.gni t.iiHp and 
direction, it can be represented by a complex number q{z), where |g(z)| is the 
fluid speed at z and Arg(9(z)) is the direction of flow at the point. 

Examples of points z and corresponding velocities q{z), for some flow, are 
shown in Figure 1.2. 


9(--l + 2i)\^ ■ 

g(l + t) = l 

q{-2 - 2i)j 
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Figure 1.2 



At a point z where the fluid is at rest (has zero speed), the flow velocity 
satisfies q{z) = 0. Such a point is called a stagnation point of the flow. 

Treating z now as a complex variable, we see that in principle the velocities at 
aU points within the fluid cross-section can be described by a complex function 
9, whose domain is the region occupied by the fluid. Before proceeding further, 
we note that the use of a complex function q to describe the fluid velocity 
pattern relies upon several significant modelling assumptions. 

First, Figure 1.1 depicts a cross-section of a flow which must in reality be 
three-dimensional, although the two-dimensional diagram is intended to show 
all the significant features of the flow. You should imagine that a ‘point moving 
with the fluid’, which at any time lies in the particular plane shown, continues 
to reside in that plane at all later times. Moreover, any cross-section of the flow 
which is parallel to the chosen plane gives a velocity pattern identical to that in 
the original plane. Such a flow is called two-dimensional. 

When we speak of a ‘point moving with the fluid’, we mean the trajectory 
taken by a smeiU volume of fluid containing the point, in the limit as the volume 
size tends to zero. It may be easiest to visualize such a ‘moving fluid point’ as 
the centre of a (very small) neutrally buoyant bead carried along by the flow. 

A second and more basic modelling assumption becomes apparent when we 
recall that a region of the complex plane is a connected open set. This implies 
that a region has no ‘local gaps’. Our model therefore assumes that the fluid 
forms a continuum, with any point of the region having a specific fluid velocity 
allocated to it. We shall assume further that the fluid velocity varies 
continuously in any direction, so that the function q is continuous. 

The third main modelling assmnption concerns what happens when the fluid 
moves on from the pattern captmred at a particular instant in Figure 1.1. In 
general, the fluid velocity q{z) at each point z will change with time. This 
complication is avoided by assuming that the flow is steady, this mesins that, 
for each fixed point z, the corresponding fluid velocity q{z) is independent of 
time. Thus Figure 1.1 represents the flow pattern at any time. 
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A ‘neutrally buoyant’ bead has 
the same density as that of the 
fluid, so that it neither rises nor 
sinks relative to the 
immediately surrounding fluid. 


Here we mean a region in the 
sense defined in Unit AS, 
Subsection 4.4. More is said 
about this below. 

The letter g is a traditional 
symbol for fluid velocity. 


According to Webster’s 
Dictionary, a ‘continuum’ is ‘an 
ideal substance or medium 
containing no vacant spaces and 
devoid of discrete structure’. In 
the context of fluids, the 
continuum hypothesis means 
that any volume, however small, 
is assumed to contain some 
matter. This assumption, and 
that of continuous velocity 
variation, describe reEil fluid 
situations accurately except at 
microscopic scales. 



Note that the assumption of steady flow is not the same as restricting the 
velocity of a point moving with the fluid to be constant as time passes. What is 
flxed is the flow velocity observed at any specified point z of the flow region. 

We now summarize the basic features of our mathematical model for fluid flows. 


Basic Mathematical Model 

We assume that 

(i) the fluid forms a continuum, and any spatial variation of the flow 
velocity is continuous; 

(ii) the flow is two-dimensional; 

(iii) the flow is steady. 

With these assumptions, we may represent the flow velocity at all times 
by a continuous complex function q, whose domain is the region occupied 
by the fluid. 


The requirements that the flow be both two-dimensional and steady may seem 
too restrictive, since there are many possible flows which are excluded. 
However, th^e are significant types of flow which satisfy these conditions, at 
least approximately; for example, the motion of air past an aeroplane wing in 
flight can for many purposes be modelled satisfactorily by a two-dimensional 
steady flow. 


Example 1.1 

Sketch the arrow pattern for the flow with velocity function 


5(z) = eW^ (^eC). 


Solution 

For each z, the flow has speed = 1 and direction given by 

Arg(e*’^/^) = 7r/4. The arrow pattern is therefore as shown in Figure 1.3. ■ 

Remark The flow in this example is called a uniform flow, because the 
velocity function is a constant function. 

Figure 1.3 was easy to draw, because all the arrows were the same length. 
When the arrow lengths (representing flow speeds) vary widely, it may not be a 
straightforward matter to sketch this type of picture. For this reason such 
diagrams are often drawn with arrows of fixed length, representing the 
direction of the flow at each point but not variations in speed. You are asked to 
adopt this approach in the next problem. 



\ 

\ 

^ / 


/ / 

/ / . 




^ / 

/ / 

/ / 


Figure 1.3 


Problem 1.1_ 

Sketch the arrow pattern, using arrows of fixed length, for each of the following 
velocity functions. 

(a) q{z) = z (zeC) 

(b) q{z) =i/z (z 6 C - {0}) 


Determining the velocity function is a main aim for any fluid flow problem 
posed within the Basic Mathematical Model. As you will see, the velocity 
function can then be used to obtain a visual representation of the actual 
motion of the fluid. 
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Here t represents time. 


To be more precise, the trajectory followed by a point moving with the fluid 
forms a path F : 'y{t) {t E I), which is called a streamline. At any point 'y{t) on 
the streamline, the velocity is g( 7 (t)) and this must be equal to 7 '(t), the rate 
of change of position, which is a tangent vector to the streamline (see Unit A4, 
page 32). Therefore we make the following definition. 


Definition a streamline through the point zq, for a flow with velocity 
function q, is a smooth path F : 7 (t) {t G I) such that 

(a) 7'(t) = qilit)), for t € /; 

(b) zo = 7(to), for some to € I- 

If q{zo) = 0 (that is, if zq is a stagnation point), then {zq} is a 
degenerate streamline, with constant paxametrization 

7(t) = 2o (f e f). 


Remarks 

1 Any such path F : 7 has infinitely many equivalent parametrizations which 
diflFer from 7 only by a time translation (that is, by having the image zq at a 
different value of to)- This corresponds to the fact that, for a steady flow, a 
description of the flow pattern is unaffected by the moment at which you 
choose to set a clock to measure time. 

2 Note that streamhnes need not be straight Unes. In the mathematical 
context, they are also called flow lines. 

As an example, consider the flow velocity function q{z) = {z G C), from 
Example 1.1. In this case. Condition (a) is 

I'it) = qiyit)) = 

so a streamUne through the point zq is given by the parametrization 
7(t) = 20 + {t G R). 

To return to the general situation, note that we may use Condition (a) to 
sketch streamlines. At each point z, the velocity q{z) is a tangent vector to a 
streamline which passes through z (unless z is a stagnation point). For 
example, we constructed the streamline pictmre in Figure 1.4 by ‘joining up the 
arrows’ of the flow pattern in Figure 1 . 1 . We indicate the flow direction along 
each streamhne by a single arrowhead. 



Similarly, the streamhnes for the flow of Example 1.1 are shown in Figure 1.5. 

Problem 1.2_ 

Sketch a few streamlines for each of the flows whose velocities are given in 
Problem 1 . 1 . 


The term ‘smooth path’ and the 
interpretation of 7 '(t) as a 
tangent vector to the path F : 7 
at the point 7 (t) were 
introduced in Unit A4, 
Subsection 4.1. 


Here we have to = 0. 



Figure 1.5 





1.2 Developing the model 

We have assumed that we can restrict attention to the steady, two-dimensional 
flow of a fluid forming a continuum, for which the fluid velocity throughout its 
flow domain is represented by a continuous complex function q. We now make 
two further modelling assumptions in order to represent the physical situation 
with a greater degree of accuracy, and to make more specific the mathematical 
problem of finding the flow velocity function q for given circumstances. 

In order to describe these restrictions, we need to introduce the ideas of 
circulation and flux. As a necessary preliminary to this, we derive a useful 
expression for the component of velocity in a given direction and introduce a 
particularly convenient type of parametrization. 

Suppose that the given direction is at an angle 6 anticlockwise from the 
positive real axis, so that it may be specified by the complex number e*®. Here 
e*® acts as a unit vector in the chosen direction, since |e‘®| = 1. 

Let q{z) be the complex number representing the velocity at a point z. 

Figure 1.6 shows two cases of how qe{z), the component of q(z) in the direction 
specified by e*®, is obtained by projecting q(z) perpendicularly onto that 
direction. 



(a) ge(z) = d>0 (b) gg(z) =-d < 0 

Figure 1.6 


We now seek an expression for the component qe{z) in terms of q{z) and 6. If 
we rotate both the complex number q{z) and the 6 direction fine clockwise 
about z through the angle 6, then q{z) goes to q{z)e~'^ and the 9 direction line 
becomes parallel to the real axis. This is shown in Figure 1.7 for the situation 
in Figure 1.6(a). Now the component of g(z)e~’® in the positive x-direction is 
Re(9(z)e-*®), 

and this is the same number as the component qe{z). Hence 
qg{z) = Re( 9 (z)e-*®). 

This last equation gives an expression for the component qg{z) of q{z) in the 
direction specified by e‘®. It tiurns out that a variation of this formula is more 
useful. This variation depends on the fact that the real parts of a complex 
number and of its conjugate are the same, so that we have 

qg{z) = Re(g(2)e-i®) = Re(q(z)e*®) . (1.1) 

Example 1.2 

Find the component of q{z) = 2i (at any point z) in the direction specified by 
eW6. 

Solution 

From Equation (1.1), with q{z) = 2i and 6 = 7r/6, we have 
q„/6{z) = Re^^e’’"/®) 

= Re(-2i(cos^+<smI)) 

= 2sin^ = l. ■ 


The component of a vector in a 
given direction is its (scalar) 
perpendicular projection onto 
that direction. 


Note that the component qe{z) 
may be positive, negative or 
zero, depending on whether the 
smallest angle between q{z) and 
e‘® has magnitude less than, 
greater thm or equal to 7r/2, 
respectively. 



_ 

d = rie(9U)c-") 
Figure 1.7 





Problem 1.3___ 

Find the component of q{z) = (at any point z) in the direction 

specified by 

Problem 1,4_ 

Use Equation (1.1) to show that the component of velocity q{z) in the direction 
at angle w/2 clockwise to that of e*® can be expressed as 


Equation (1.1) and the result of Problem 1.4 can be rewritten in terms of the 
conjugate velocity function g, which has the same domain as q and rule 

q{z) = q{z). 

We have 

qe{z) = Re(g( 2 :)e*®) and q(e-^/2){z) = Im(g(z)e’®), (1.2) 

which will be of use shortly. 

Now suppose that F : 7 (t) (t € /) is a smooth path of finite length lying in the 
domain of a flow velocity fimction q. We wish to define two quantities, 
circulation and fiux, associated with F and q, and it will be convenient to 
assume that 7 is a unit-speed parametrization; that is, 

l 7 '(t)| = 1 , for t € /. 

The following result shows that we can always do this. 


Theorem 1.1 if F : 7 ^( 4 ) {t e [a, 6]) is a smooth path of length L, then 
there is an equivalent parametrization 7 ( 3 ) (s e [ 0 ,L]) such that 

| 7 '(s)| = 1, for 0 < s < L. (1.3) 


Proof Since 7 ^ is smooth, the real function 1 1 —* | 7 l(i)| is continuous on [a, 6 ], 
and so it has a primitive h, say. Assuming, as we may, that h{a) = 0, we obtain 

h{b) = h{b) - h{a) = j \‘y[{u)\du = L, 

by the Fundamental Theorem of Calculus. Because h'{t) = | 7 i(t)| > 0, we 
deduce that the function 7 ( 3 ) = 7 i(h“^(s)) (s e [ 0 , X]) is an equivalent 
parametrization with 

| 7 '(s)| = |7i(*~^(®))(^~^)^(^)| (Chain Rule) 

= \'yi{t)\/h'{t) (Inverse Function Rule : t = /i~^(s)) 

= 1 , for 0 < s < L, 
as required. ■ 

Rsmark if Equation (1.3) holds, then the parameter s measures length along 
F since, if 0 < sj < S 2 < X, then the length of F from 7 (si) to 7 ( 32 ) is 

|7'(s)|ds = ^ Ids = 32 - 31 . 

Note that it is traditional to use the letter s as parameter with unit-speed 
parametrizations. 


For example, if q(z) = iz, then 
q{z) = = -iz. 


Note that F is not necessarily 
part of a streamline. 


For example, 

7(t) = 2(cosit-hisinit) 

(t € [0,47r]) 

is a unit-speed parametrization 
of the circle jzj = 2. 


Unit Bl, Frame 5 


_ 1 
- h'{ty 
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The use of unit-speed parametrizations makes it possible to define the 
tangential and normal components of a flow velocity function g at a point on a 
smooth path F. 


Definitions If r : 7(3) (s € [0, L]) is a smooth path with unit-speed 
paxametrization, which Ues in the region 7?. of a flow velocity function q, 
then, for each s G [0, L], the flow velocity 9(7(3)) has 

(a) tangential component qris) in the direction specified by 7 ^( 3 ); 

(b) normal component gjv(s) in the direction specified by 


These components axe shown in Figure 1.8 in a case when both axe positive. 



Note that 7'(s) is a unit 
tangent vector at the point 7(3) 
on r, for each s € [0, L], and 
—i7'(s) is a unit normal vector. 


We now define circulation and flux in terms of these velocity components. 


Dsfinitions Let r be a smooth path of length L in the domain of a fiow 
velocity function q, and let 7 : [0, L] —> F be the equivalent unit-speed 
paxametrization of F. 

(a) The circulation of q along F is 


Cr= qT{s)ds. 

Jo 

(1.4) 

(b) The flux of 9 across F is 




•^r = / qN{s)ds. 

Jo 

(1.5) 


Note that these are real 
integrals, so both the circulation 
Cr and the flux Tr are real 
quantities. 


In words, the circulation of q along F is 
L X the average value of qx on F, 
and the flux of q across F is 

L X the average value of q^f on F. 

(You will be asked to calculate a circulation and a flux in Problem 1.7.) 

Now for each s € [0, L], 7 ^( 3 ) has modulus 1, and so we have 7 ^( 3 ) = e*® for 
some 6. Thus 


In physical terms, flux is a 
measure of the overall rate at 
which fluid is crossing F. 


Qt{s) = 93 ( 7 ( 5 )) and 
so that, by Equations ( 1 . 2 ), 
9t(s) = Re(9(7(3))7'(3)) 


9JV(S) = g(8-7r/2)(7(5)), 

and q^is) = Im( 9 ( 7 (s)) 7 '(s)). 


( 1 . 6 ) 
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It follows from the definitions of circulation and flux, in Equations (1.4) 
and (1.5), that 

Cr + iJ^v = [ (9t(s) + *9jv(s ))ds 
Jo 

^ Jo 

that is, 

Cr + i^T = J q{z)dz. (1-7) 

In deriving this formula, we assumed for simphcity that the path F has a 
unit-speed parametrization. However the formula is valid also with any 
equivalent parametrization for the path. Furthermore, this result for smooth 
paths can be extended easily to a corresponding result for contours. Indeed, 
Equation (1.7) holds for any contour F, provided that the circulation Cr is 
defined as the sum of the circulations along the smooth paths which form F, 
and the flux Tr is defined as the sum of the fluxes across these smooth paths. 

We are now almost ready to make the final modelling assumptions. These are 
based on the following definitions. 


Definitions a velocity function q with domain a region TZ is 

(a) locally circulation-free if Cr = 0 for each simple-closed contour F 
in TZ whose inside also lies in TZ; 

(b) locally flux-free if .Fr = 0 for each simple-closed contour F in 71 
whose inside also lies in TZ. 


Remarks 

1 Figure 1.9 shows such a contour F in a region TZ which is not 

simply-connected. Note that if TZ is simply-connected, then the inside of any 
simple-closed contour in TZ automatically hes in TZ. 

2 The word ‘locally’ occurs here because, for each point z gTZ, the condition 
that the inside of the contour F should he within TZ is satisfied automatically if 
F lies in any open disc in TZ with centre at z. 

We shall assume for our model that any flow to be considered has a velocity 
function which satisfies both of the conditions just defined. 


Definition a model flow of a fluid is described by a continuous 
complex velocity function (whose domain is a region) which is locally 
circulation-free and locally flux-free. 


The locally flux-free condition is equivalent (for a constant-density fluid) to the 
Principle of Conservation of Mass, which states that mass (matter) is neither 
created nor destroyed. In the context of fluid flows this means the following: 
for any simple-closed contour F whose inside is in the flow region TZ, and at Einy 
time, any flow of fluid mass (per unit depth) inwards across F must be exactly 
balanced by a mass flow outwards, so that the mass of fluid associated with the 
inside of F remains constant. 

The need for the locally circulation-free condition is less obvious, and indeed 
there are actual flows, satisfying approximately all of the other specified 
conditions, which are not locally circulation-free. Roughly speaking, this 
condition permits us to restrict attention from all the possible flows for a given 
situation to the unique flow among them which has ‘minimum energy’ 
associated with it. In fact, in terms of a small neutrally buoyant bead, whose 
surface velocity matches that of the adjacent fluid, the locally circulation-free 
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This is a consequence of 
Unit Bl, Theorem 2.1. 



Figure 1.9 


In vector c^llculus notation, the 
locally flux-free condition is 
expressed as 
div q = 0. 


In vector calculus notation, the 
locally circulation-free condition 
is expressed zs 
curl q = 0. 

By a ‘given situation’ here, we 
mean the specified configuration 
of physiceJ boundaries for the 
flow, limiting flow speed at large 
distances, and so on. 




condition means that the bead will not rotate. In other words, the flow has no 
‘swirl’ or ‘local angular velocity’. For this reason, a flow with this property is 
also called irrotational. 

It follows immediately from Equation (1.7) that 

a fluid flow is a model flow on a region TZ if and only if its velocity 
function q satisfies the condition 

j^q{z)dz = 0, 

for each simple-closed contour F in 7?. whose inside also Ues in TZ. 

This important fact will be apphed in the next subsection. 


1.3 The conjugate velocity function 

Our model flow is described by a continuous velocity function q, whose domain 
is a region TZ. As pointed out in the previous subsection, the conjugate velocity 
function q may be defined by 

q{z) = q{z) (zeTZ). 

The function q is continuous, since q is continuous. 

Now, by oiu: assumptions, q is locally both circulation-free and flux-free. 

Hence, as you saw above, we have 

j^q{z)dz = Q, (1.8) 

for each simple-closed contour T in 72. whose inside hes in 72. In particular. 
Equation (1.8) holds for all rectangular contours whose insides he in 72, and so 
it foUows from Morera’s Theorem that the conjugate velocity function q is 
analytic on 72. 

The converse of this result also holds. If g is any function einalytic on a region 
72 then, by Cauchy’s Theorem, 

J^g{z)dz = 0, 

for each simple-closed contour T whose inside lies in 72. By Equation (1.7), the 
function q defined by q{z) = g{z) must then be locally both circulation-free and 
flux-free. We have thus estabhshed the foUowing result. 


Theorem 1.2 A steady two-dimensional fluid flow with continuous 
velocity function g on a region 72 is a model flow if and only if its 
conjugate velocity function q is analytic on 72. 


This theorem provides a method for finding a profusion of model fluid flows, 
since the conjugate of each analytic function provides a corresponding model 
flow velocity function. 

Example 1.3 

Find the model flow velocity function corresponding to the cinalytic function 
giz) = {z e C). 

Solution 

The conjugate velocity function is q{z) = so that the model flow 

velocity function is 

q{z) = 6"/“ {z e C). ■ 


In vector calculus notation, the 
condition for the imaginary part 
of this equation is equivalent to 
div q = 0, via Gauss’ Divergence 
Theorem. The condition for the 
real part is equivalent to 
curl q = 0, via Stokes’ Theorem 
(or Green’s Theorem in the 
Plane). 


Note that 


Unit B2, Theorem 5.4 


Unit B2, Theorem 1.2 


The streamlines for this flow are 
shown in Figure 1.5. 
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Problem 1.5___ 

For each of the following analytic functions g, find the corresponding model flow 
velocity function. Draw a fixed length arrow pattern to illustrate each answer, 

(a) g{z) = z (2 G C) 

{h)g{z) = l/z (^gC-{0 }) 

{c) g{z)^i/z (^eC-{0 }) 

Problem 1.6___ 

Which of the velocity functions specified in Problem 1.1 represents a model 
fluid flow? 


Problem 1.7_ 

Let r be the unit circle {z: \z\ = 1}. 

(a) Use Equation (1.7) to find the flux across T of the velocity function which 
you found in Problem 1.5(b). 

(b) Use Equation (1.7) to find the circulation along F of the velocity function 
which you found in Problem 1.5(c). 

(c) Explain why your answers to parts (a) and (b) do not contradict 
Theorem 1.2. 


Problem 1.7 demonstrates that the flux or circulation of a flow may be 
non-zero for certain closed contours even when the flow is locally both 
circulation-free and flux-free. The example of Problem 1.5(b) features a source 
of strength Tj- = 27r at the origin, corresponding to outward flow, whereas that 
of Problem 1.5(c) features a vortex of strength Cr = -27r at the origin, 
corresponding to clockwise flow. More generally, we have the following 
definitions, which are illustrated by Figure 1.10. 


DofinitionS Let g be a model flow velocity function with domain a 
region K, and let D be a punctured open disc in "R with centre a. Then 

(a) a is a source of strength Tii !Ft = T ^0 lox each simple-closed 
contour F in D which surrounds a; 

(b) a is a vortex of strength C if Cp = C ^ 0 for each simple-closed 
contour F in £) which surrounds a. 


Remarks 

1 If the point a is a source or a vortex, then it is not in the region R. In fact, 
q then has an isolated singularity at a {Unit B4, Subsection 1.1). 

2 The fact that .Fp and Cp are independent of the choice of the simple-closed 
contour F is a consequence of the Shrinking Contour Theorem (Unit B2, 
Theorem 1.4). 

3 In these definitions, F and C may be positive or negative. A source of 
negative strength is called a sink, since fluid is then flowing in to the point a 
rather than out from it. Negative vortex strength indicates clockwise rather 
than anticlockwise flow. 

To conclude this section, we ask you to use the Cauchy-Riemann equations to 
obtain an alternative formulation of the condition that a model flow velocity 
function should be locally both circulation-free and flux-free. 



Figure 1.10 


Unit A4, Subsection 2.1 
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Problem 1.8- 

(a) Let q be a continuous complex velocity function with domain a region H, 
and suppose that q\ = Re g and 92 = Im g have partial derivatives with 
respect to x and y which are continuous on TZ. 


Prove that g is a model flow velocity function on TZ if and only if 

^ + ^=0 and ^-$1 = 0 onK. 
dx dy dx dy 

(b) Use your result from part (a) to verify that the function 


q{z) = i/z (z e C - {0}) 


of Problem 1.1(b) represents a model flow. 


Here we eire treating gi and q2 
as functions of (x, y) rather than 
of z = X + fy. 


Remark The two equations in Problem 1.8(a) can be used to show that if g is 
a model flow velocity function, then both the following equations hold: 

9x2 + Qy2 0 Q^2 ^ ay2 

Thus both gi and g 2 satisfy the second-order partial differential equation 
d'^u d'^u 


- = 0 on 7?.. 


9x2 gy2 


= 0 , 


Here d^qi/dx^, for example, is 
the second parti2il derivative of 
gi with respect to x, that is, the 
second derivative of gi with 
respect to x while keeping y 
constant. 


which is called Laplace’s Equation. Functions which are solutions of 
Laplace’s Equation are called harmonic. Thus the real and imaginary parts of 
any model flow velocity function are harmonic, as are the real and imaginary 
parts of any analytic function. 


Pierre Simon Laplace 
(1749-1827) was a leading 
I^ench mathematician and 
astronomer. 


You showed in Problem 1.8(a) that a function g is a model flow velocity 
function (that is, locally both circulation-free and flux-free) if and only if the 
two partial differential equations given in the problem both hold. In fact, this 
result can be spht into two constituent parts, though we shall not prove this. 


Theorem 1.3 Let g be a continuous complex velocity function on a 
region TZ, and suppose that gi = Re g and q^ = Im g have partial 
derivatives with respect to x and y which are continuous on TZ. Then g is 
(a) locally circulation-free on H if and only if 


9q2 dqi ^ 
dx dy 


on TZ-, 


(b) locally flux-free on 7^ if and only if 


9qi d^ _ 
dx dy 


on TZ. 


In vector calculus notation, 
these conditions may be written 
respectively as 

(a) curl q = 0; 

(b) div q = 0. 


These conditions on gi and g 2 are often easier to apply t)ian the original 
deflnitions of locally flux-free and locally circulation-free, especially when only 
one of these properties holds. 


Problem 1.9_ 

Use Theorem 1.3 to show that 

(a) the flow velocity function g(z) = z of Problem 1.1(a) is locally 
circulation-free, but is not locally flux-free; 

(b) the flow velocity function g(z) = iz is locally flux-free, but is not loceiUy 
circulation-free. 


In the next section you will see further consequences of the close link between 
model fluid flows and analytic functions. 
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2 COMPLEX POTENTIAL FUNCTIONS 

After working through this section, you should be able to; 

(a) understand and apply the definition of the complex potential function', 

(b) find families of streamlines, using the complex potential function or stream 
function; 

(c) recognize complex potentials for a source, sink or vortex fiow; 

(d) appreciate the main features of a flow involving a doublet, and the effect of 
placing a doublet in a uniform flow. 


2.1 An equation for streamlines 

For the rest of this unit we deal only with model fluid flows, for which 
Theorem 1.2 provides a fruitful connection to the theory of analytic functions. 
You have already seen how this link can be applied either to generate model 
flow velocity functions or to test, using the Cauchy-Riemann equations, 
whether a velocity function does indeed describe a model fluid flow. 

We look next at the application to fluid flows of another important result for 
analytic functions, namely, the Primitive (or Antiderivative) Theorem, which 
states that if a function is analytic on a simply-connected region TZ, then it has 
a primitive on TZ. 

Suppose that g is a model flow velocity function with domain TZ. Then, by 
Theorem 1.2, the conjugate velocity function q is an analytic function, with the 
same domain TZ. Now TZ may not be simply-connected, but if 5 is a 
simply-connected subregion of TZ then, according to the Primitive Theorem, q 
has a primitive function Q, say, on 5. In other words, there exists a function fl, 
analytic on its domain S, such that 

= q{z), for 2 e 5. 

In the context of model fluid flows, a function Q which is a primitive of q is 
called a complex potential function for the flow. By taking the complex 
conjugate of the equation above, we see that the flow velocity function q is 
given in terms of the complex potential function by 

q{z) = Cl'{z), for 2 6 5. 

The significance of the complex potential function requires some explanation. 
First, note that, from Equation (1.7), we have 

Cr + iTr = j q[z) dz, 

where Cr and .Fr are, respectively, the circulation along and flux across any 
contour F in TZ. Now, provided that F lies within the simply-connected 
subregion 5 on which the complex potential Q is defined, we may apply the 
Fundamental Theorem of Calculus to deduce that 

Cr + iTr = J ^'(2) dz = il{P) - fl{a), (2.1) 

where a and P are the initial and final points of F. The real and imaginar y 
parts of this equation are, respectively, 

Cr = Re(n(/3)) - Re(n(Q)) and JTr = Im(n(/?)) - Im(n(a)). (2.2) 

These formulas show that, for a model flow velocity function, the circulation 
along (or flux across) any contour within 5 from a to /3 is the same. (This is an 
example of the Contour Independence Theorem {Unit B2, Theorem 1.3).) For 
example, for the contours shown in Figmre 2.1, the integral of f2'(2) along the 
contour Fi and that along F2 are both equal to fi(/I) - ^(q). 


Unit B2, Theorem 5.3 


Often this is abbreviated to 
‘complex potential’. The use of 
the Greek letter O (capital 
omega) for the complex 
potential function is traditional 
in the study of fluid dynamics. 


Note that F need not be a 
closed contour here. 

Unit Bl, Theorem 3.1 



Figure 2.1 
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Now flux was deflned (Equation (1.5)) as the integral along a smooth path F : 7 
of the normal component 9 n(s) of the flow velocity qi'ris)). Thus if the flow 
velocity at each point 7(5) of F is a tangent vector to the path (so that 
qff{s) = 0), as shown in Figure 2.2, then there is zero flux across F. Paths of 
this type are the trajectories of points moving with the fluid. They were given 
the name ‘streamlines’ in Subsection 1.1, where we indicated their usefulness as 
a means of viewing a fluid flow. 



7(53) 9(7(83)) 


Figure 2.2 


A streamline F therefore has the property that there is no flux across it. Thus 
if a and /3 both lie on a streamline F, then, by applying Equation (2.2) to the 
part of F which joins a to 0, we have 

Im(n(/?)) = Im(n(a)). 

Thus we have proved the following result. 


ThGOrorn 2.1 if n is a complex potential function on a simply-connected 
subregion S of the flow region TZ, then 

Im(n(2r)) is constant along each streamline within S. 


Remarks 

1 This theorem states that, for each streamUne F in 5, there is a constant k 
such that Im(n(jz)) = fc for all z € F. However, the set {z e 5: Im(n(z)) = k} 
may comprise more than one streamhne. You will see examples of this type in 
Subsection 2.2. 

2 The result refers to the portion of a streamline F within a simply-connected 
subregion S of the flow region TZ. If F extends beyond <S, then it is possible to 
cover the extension by an appropriate sequence of simply-connected subregions 
Sn with associated complex potentials n„ (n = 1,2,...), where Si = S and 

Hi = 0 (see Figure 2.3). Each complex potential of the sequence after the first 
is related to the previous one by direct analytic continuation, and the constant 
value of Im on F will then be the same for each n. 

Iin(n.(z)) = k 



3 Theorem 2.1 can also be used to show that each point of the flow region has 
just one streamline through it. Despite this, it may sometimes appear that 
there are two streamlines passing through a point. This impression can occur 
only at a stagnation point, which is itself a degenerate streamline. 


Unit CS, Subsection 1.1 
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Theorem 2.1 indicates one reason why complex potential functions are of 
interest. If 


^l{z) = $(x, y) + i^{x, y), where z = x + iy, 

is an analytic function on a simply-connected region 5, where ’i' are 
real-valued functions, then the family of cmves described by the equation 

'if{x,y) = constant 


The use of the Greek letters $ 
and ^ (capital phi and psi) for 
the real and imaginary parts of 
n is standard practice in fluid 
dynamics. 


forms the streamlines of a model fluid flow. The function 4' = Im f2 is therefore 
called the stream function for the flow. 


In the next subsection we look at several model flows, given by flow velocity 
functions. In each case we find a complex potential function and then the 
corresponding stream function. 

Problem 2.1___ 

Show that the flow velocity function q may be expressed directly in terms of 
the stream function ^ as 




This gives the flow velocity q{z) 
at z in terms of the stream 
function 


Remark For a non-model flow whose velocity function is locally flux-free (but 
not locally circulation-free), there still exists a stream function ’4', defined in 
terms of the velocity function g as in the result of Problem 2.1, for which the 
streamlines are given by ’^'(x, y) = constant. 

However, in this case there is no anal3rtic conjugate velocity function and hence 
no complex potential; that is, there is no analytic function fl for which 
^ = Im$l. Such flows are not considered further in this unit. 


2.2 Examples of fluid flows (audio-tape) 

This audio tape presents several examples of model fluid flows, some of which 
you have seen earlier. In each case the starting point is a model flow velocity 
function, from which we derive a complex potential and the corresponding 
stream function. This enables the streamline diagram to be drawn. Despite the 
simplicity of these cases, you will see later that they can be used as the 
building blocks for more comphcated model flows. 

NOW START THE TAPE. 
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Problem 2.3- 

(a) Let h>0. Show that the model flow velocity function 



has complex potential function 

Cl{z) = Logz - Log(2: — h) (2 e C — {x e IR : X < h}). 

This represents the flow resulting from the combination of a source at the 
origin and a sink a.t z = h, both of strength 27r. 

(b) Use elementary geometry to show that the streamlines of this flow axe 
circular axes terminating at z = 0 and z = h, with the centres of the circles 
lying on the line x= \h. 


Problem 2.4_ 

The model flow velocity function 


represents the same type of flow as in Problem 2.3, but with the strengths of 
source and sink now equal to 2-K/h. Show that, if the hmit as /i —♦ 0 is taken 
(in which process the strengths of source and sink axe always inversely 
proportional to the distance between them), then the resulting flow is that 
arising from a doublet. 



3 FLOW PAST AN OBSTACLE 

After working through this section, you should be able to: 

(a) explain the main features of the flow past a disc due to a uniform stream 
with circulation; 

(b) explain what is meant by an obstacle, and state the Obstacle Problem; 

(c) solve the Obstacle Problem in appropriate cases by applying the Flow 
Mapping Theorem, for flows with or without circulation around the 
obstacle; 

(d) explain what is meant by a Joukowski aerofoil. 

3.1 Flow past a circular cylinder 

In this subsection we return to the flow considered in Frame 7 for which the 

model flow velocity function is 

q{z) = l-a^/z^ (^6C-{0}), (3.1) 

and a corresponding complex potential function is 
n{z)^z + a^/z. 

Here, and throughout this section, the symbol a denotes a positive real number. 

Problem 3.1_ 

(a) Verify that the points z = a and z = —a are the only stagnation points of 
the flow deflned by Equation (3.1). 

(b) Show that the flow defined by Equation (3.1) satisfies 

q{ae’^) — (-2isint)e*‘, for 0 < t < 27r. 

Hence derive the flow speed at the point z — ae‘‘, on the circle {z : \z\ — a}. 


In the audio tape (Frame 8) we claimed that Equation (3.1) can be used to 
describe the flow due to a uniform stream past a solid cylinder of circular 
cross-section (see Figure 3.1); we shall denote this cross-section as follows: 

Ka = {z: 1^1 < o}, where a > 0. (3.2) 

This claim was based on the fact that dKa — {z : |2:| = a} is made up of 
streamlines. However, there are other possible model flows for which dKa is 
made up of streamlines. The complex potential function 

f2(a) =-icLogz (z e C,), 

where c e IR, corresponds to vortex flow (see Frame 4), and its streamlines are 
all circles with centre 0. Thus it seems reasonable to suppose that the complex 
potential function 


^a,c{z) = Z-\ -icLogZ (z € Cfl.) 


(3.3) 


should correspond to a model flow velocity function for which dKa is made up 
of streamlines. This flow will play a major role in this unit and so we record its 
main properties in a theorem. 


Indeed, 

WiT) = q{z). 


Here the parameter t does not 
stand for time, as it did in 
Section 1. 



Figure 3.1 


This complex potential is c 
times that for the vortex flow 
Frame 4. 
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Theorem 3.1 For a > O, c € R, the model flow velocity function 

= (^eC-{0}), (3.4) 

corresponding to the complex potential function f2o,c of Equation (3.3), 
satisfies the following properties: 

(a) lim qa,c{z) = 1; 

(b) dKa is made up of streamlines for qa,c\ 

(c) for any simple-closed contour F surrounding Ka, 

(i) Cr = Re J q^{z) dz = 27rc, 
and 

(ii) .fr = Im j '^(z) dz = 0. 


Remarks 

1 The expression lim qa,c{^) is defined in a similar way to an ordinary limit 
as z tends to a, where q € C; in particular, hm qa,c{z) = 1 means that 

if {zn} is a sequence such that z„ —> oo, then qa,c{zn) —* 1- 

2 Since Cr = 27rc for any simple-closed contour F surrounding Ka, we call 27rc 
the circulation around the disc Ka- 

3 In fact it can be shown that Property (c)(ii) is a consequence of Property (b). 


Proof of Theorem 3.1 


Property (a) follows immediately from Equation (3.4). To verify Property (b), 
note that, for z € dKa, and z -a, we have 

Im(fia,c('2)) = Ini(2 -(- a^/z) -I- Im( —icLogz) 

= -clogg a (Frames 4 and 7). 

Since —clogg a is a constant, the set dKa, is made up of streamlines for qa,c- 
To verify Property (c), we use the Residue Theorem to obtain 


/r’“•■(") * = /r(‘-?-7)* 
= 27ri(-ic) = 27rc. 


Hence Cp = 27rc and =0. ■ 


The additional vortex term (which gives rise to the non-zero circulation 
around F) affects the positions of the stagnation points, as we now ask you to 
verify. 

Problem 3.2_ 

It follows from Equation (3.4) that qa,c{z) = 0 if and only if 
z^ — icz — =0. 

By solving this equation, find the stagnation points of the flow. Take c to be 
negative and describe the locations of these stagnation points when 
(a) -2a < c < 0; (b) c = -2a; (c) c < -2a. 


Theorem 3.1 allows us to think of qa,c as the model flow velocity function 
obtained by placing the ‘obstacle’ Ka into a uniform stream (for which 
q{z) = 1) and requiring that there is no flow of fluid across any segment of 
dKa, but that the circulation of the flow around Ka is 27rc. 


Note that 

fCW=9a.c(^) 
on C,r, but qa,c has domain 

C-{0}. 


Unit AS, Section 3 


To include z = —a, we consider 
a logarithm function with a 
different cut, as discussed on the 
audio tape. 


The stagnation points for the 
case c = 0 were found in 
Problem 3.1. 


We ask you to concentrate on 
the cases with negative c 
because these will be 
particularly relevant later in the 
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Each streamline of this flow satisfies an equation of the form 
Im(na,c(z)) = Im(z + 0^ jz) + Im( —ic Log z) 

= constant. 

Unfortunately it is not possible to solve these equations (for y in terms of x, or 
vice versa) when c / 0. The streamline diagrams (for various values of c < 0) 
in Figure 3.2 are based on computer plots. The directions of the arrows can be 
obtained from Equation (3.4); a quick method is to start with q{z) = 1 when 
|z| is large, and then use continuity. 




Figure 3.2(d) should be familiar 
— it is the cover design of 
Block D! 


Notice how the shapes of these streamlines are closely related to the locations 
of the stagnation points in each case. For c = 0 (Figure 3.2(a)), the two 
stagnation points are symmetrically placed on the real axis at z = ±a, by 
Problem 3.1. As c decreases from 0 to -2a, the two stagnation points 
(obtained in Problem 3.2) 

z = ^ic ± Y a^ — ^c^ 

move round the circle |z| = a, remaining symmetric with respect to the 
imaginary axis (Figure 3.2(b)), until c = —2a, when they coalesce to form a 
single stagnation point at z = -ia (Figure 3.2(c)). For c < -2a, the single 
stagnation point moves off the circle down the imaginary axis (Figure 3.2(d)) at 


The other stagnation point lies 
inside the open disc 
{z : |z| < a}, as shown in the 
solution to Problem 3.2(c). 


. = i(ic-v'i?3^), 

A similar analysis can be carried out for positive values of c, resulting in 
streamline diagrams which are reflections in the real axis of the diagrams in 
Figure 3.2. 
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3.2 The Obstacle Problem 

We now discuss the problem of determining the flow past a more general 
obstacle in a uniform stream. To begin with, we make precise what is meant 
here by an obstacle: this should be thought of as the two-dimensional 
cross-section of an object, extending indeflnitely fax perpendicular to the 
complex plane. 


Definition An obstacle is a compact connected set K in C, such that 
C — if is also connected. 


Remark it follows that C - if is a region. 

According to this deflnition, the closed disc Ka = {z: \z\ < a} is an obstacle, as 
is the union of an elhpse with its inside, blit the open disc {z : \z\ < a} is not 
an obstacle. Also the line segment joining a and /3 is an obsteicle. Later in the 
unit this hne segment is denoted by L{a,(3); thus 

Lia,P) = {{l-t)a + tl3:0<t<l}. 

Note that L{a,0) and L{l3,a) describe the same line segment. Also, if a, /3 are 
both real, and a < 0, then L{a,0) = L{P,a) = [a,0\. 

The Obstacle Problem is to And, for a given obstacle if and a given c € R, a 
model flow velocity function which, loosely speaking, satisfles Properties (a), 

(b) and (c) of Theorem 3.1. Properties (a) and (c) are easy to state for a 
general obstacle, but Property (b) is more tricky. We want dK, the boundary 
of if, to look, as far as possible, as if it is made up of streamhnes, but it is 
diflScult to express this in terms of the velocity function q, because q may 
behave badly near some points, such as comers of dK. Instead we express this 
streamUne condition in terms of a stream function, which is the imaginary part 
of a complex potential function fl for q. As with the complex potential function 
fla.c) used in the proof of Theorem 3.1, we cannot expect that fl will be 
cinalytic on the whole of C — if, but we try to make its region of analyticity as 
large as possible, so that we can give meaning to the streamline condition 

Im(n(z)) = constant, for z £ dK, 
even though dK may not be in the domain of and q. 


Obstacle Problem 

Given an obstacle K and a real number c, we seek a model flow velocity 
function q with domain the region 'R = C — K, satisfying the following 
properties: 

(a) ^1^9(2:) = 1; 

(b) there is a complex potential function f2 for q on either K ox Tt - T,, 
where E is a simple smooth path in 'll, and a real constant k such that 

for e£ich a £ dK, we have lim Im(f2(2)) = fc; 

(c) for any simple-closed contour P surrounding if, 

Cr = 27rc. 


The notions ‘compact’ and 
‘connected’ were defined in 
Unit AS, Subsections 5.1 
and 4.3, respectively. 


In particular, we shall consider 
the flow past line segments of 
the form L{—a, a), whose centre 
is at 0. 


fla.c is analytic on 

(C - ifa) - {x e R : X < -a}, 

which is a subset of C^. 


A simple path is one which 
does not intersect itself (see 
Unit B2, Subsection 1.1, where 
the dom^un of the 
parametrization of the path was 
assumed to be of the form [a, 6]). 
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Remarks 

1 Property (b) is illustrated in Figure 3.3 for the Ccise where there is a complex 
potential for q on 7^ — E. 

Figure 3.3 



Recall that lim Im(n(2:)) = k means that a is a limit point of 7^ - E and, for 
each sequence {2„} in 7^ - E such that Zn —* a, we have Im.{fl(z„)) —► k. Thus, 
in this case. Property (b) could be expressed, equivalently, by saying that the 
function Im f2 has a continuous extension from 7^ — E to (7^ - E) U dK with 

lmn{z) = k, hrzedK. 

2 It is clear from the proof of Theorem 3.1 that the function qa,c in 
Equation (3.4) solves the Obstacle Problem for K = Ka = {z ■.\z\ < a}, with 
circulation 27rc around K. The simple smooth path E in Property (b) is the 
part of the negative real axis which lies in C - ATa (see Figure 3.4), and the 
constant k is -clogg a, as in the proof of Theorem 3.1. Note that if c = 0, then 
the simple smooth path E is not needed since the complex potential function 
0(2) = z + a^/z is analytic on C — {0} and hence on C — ifa- 

3 Despite the comments before the statement of the Obstacle Problem, we 
shall refer to flow velocities on dK when it makes sense to do so. This will be 
the case when q has a continuous extension from C — A" to dK (as does qa,c 
when K — Ka)- When this notion is ambiguous (for example, when the 
obstacle is a line segment with different limiting velocities as it is approached 
from either side) we find it convenient to refer to velocities ‘on either side of 
the boundary’. We consider such an example in Remark 2 after Example 3.1. 

4 You will have noticed that Property (c) does not include the statement 

= 0. In fact, it is a consequence of Property (b) that Jr = 0, so to add this 
to the requirements for q in the Obstacle Problem would be to introduce 
redundancy. In Theorem 3.1, the property .Fp = 0 is a consequence of the form 
of 9a,c given in Equation (3.4), even though it can be deduced from 
Property (b). 

5 The quantity 27rc in Property (c) is the circulation around the 
obstacle K. 

In the next problem we ask you to verify that the velocity function ^(z) = 1 
solves the Obstacle Problem in a simple case with zero circulation (so that the 
simple smooth path E is not needed). 

Problem 3.3_ 

Show that the model flow velocity function 
q{z) = 1 

solves the Obstacle Problem for any line segment K parallel to the real axis, 
with zero circulation around K. 


Unit AS, Section 3 4 


n-7: 



Figure 3.4 
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3.3 Flow Mapping Theorem 

Problem 3.3 shows that the Obstacle Problem has a very simple solution for an 
obstacle K consisting of a line segment parallel to the real axis when the 
circulation around K is 0. However, if the line segment is not parallel to the 
real axis, or if the circulation is not 0, the Obstacle Problem becomes much 
harder to solve. 

The velocity function q{z) = 1, which solved the Obstacle Problem in 
Problem 3.3, is one which you might have guessed. Rather than attempt to 
guess solutions to the Obstacle Problem for particular obstacles K and given 
circulations, we shall use a technique based on conformal mappings. The idea 
is that if X is a given obstacle and / is a suitable conformal mapping from 
C — if onto C — ifa, where Ka — {z'- 1^1 < a}, then it is possible to use our 
solution qa,c of the Obstacle Problem for Ka in order to solve the Obstacle 
Problem for K. The details are given in the following result. 


Theorem 3.2 Flow Mapping Theorem 

Let K be an obstacle and let / be a one-one conformal mapping from 
C — if onto C — ifa, where a > 0, such that 

f{z) = z + ao + ^ + ^ + ---, for I 2 I > R, (3.5) 

where R> 0 and ao,a_i,a_ 2 ,... € C. Then the velocity function 

q{z) = 9a.c(/(^))7^ (^ € C - if) (3.6) 

is the unique solution to the Obstacle Problem for if, with complex 
potential function 

n = Ha.c o /. 


The proof of this theorem is 
given at the end of this 
subsection. 


ifa, 5o,c and fJa.c are given in 
Equations (3.2), (3.4) and (3.3), 
respectively. 


Remarks 

1 Theorem 3.2 is illustrated in Figure 3.5, which also indicates that the 
conformal mapping / maps streamlines to streamlines. This is so because 
^ o / and so, for any constant fc, the image under f of the set 

{z : Im(H( 2 )) = k} = {z: Im(n„,,(/( 2 ))) = fc} 
is the set {lu : Im(na,c(t'^)) = k}. 

Similarly the conformal mapping maps streamlines to streamlines. 



Figure 3.5 


2 It follows from Equation (3.5), which gives the form of the Laurent series 
about 0 for / on {2 : | 2 | > R}, that for large | 2 |, / behaves approximately like a 
translation. Hence, / cannot be a non-trivial rotation or a scaling. 

Theorem 3.2 reduces the Obstacle Problem for K to the problem of finding a 
one-one conformal mapping from C - ff onto C — Ka, for some a > 0. In 
general, this is a hard problem, but there are some cases which can be solved 
explicitly. One simple example is the subject of Problem 3.4. 
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Problem 3.4___ 

Use the Flow Mapping Theorem to solve the Obstacle Problem for the disc 
K = {z-.\z-^\<a}, 

where (3 is any complex number and the circulation around the disc is 27rc. 


Another example employs the inverse function of the Joukowski function. In 
Unit Dl, Subsection 4.4, you saw that the Joukowski function 

J{z) = z + ^ {z£C- {0}) 

gives a one-one conformal mapping from C- Ki, where Ki = {z: \z\ < 1}, 
onto C- K, where K is the line segment L(-2,2) = [-2,2]. Therefore 
which is the inverse function of the restriction of J to C - ifi, gives a one-one 
conformal mapping from C — K onto C - Ki (see Figure 3.6). The rule for J~^ 
was found to be 

= 5(w + u;v'1-4/u;2) . 

Also, since 

dKi == {e*‘ : 0 < t < 27r} , 
we have 

J{dKi) = -I- l/e*‘: 0 < t < 2?:} 

= {2 cos t: 0 < t < 27r} = [-2,2]. 



w-plane 




Figure 3.6 


It follows that, if we reverse the roles of z and w here, then we can use as 
the function / in the Flow Mapping Theorem and so solve the Obstacle 
Problem for the line segment K = [-2,2] with non-zero circulation. 


Example 3.1 

Use the function J~^ to show that the solution to the Obstacle Problem for 
K = [—2,2] with circulation 27rc around K is 

"-(it*!?) 

Solution 


As pointed out above, 

f{z) = J~'^{z) = I (z + zy/l - 4/z^^ 


is a one-one conformal mapping from C — K onto C — Ki, so that a = 1 in this 
case. To verify Condition (3.5), we need to express J~^ as a Laurent series 
about 0, valid on the exterior of some closed disc. To do this we use the 
binomial series 


(1 - 4/z^f^ = 1 + l(-4/.^) +(1)(-1) + . 


for 1^1 >2, 


for |—4/z^l < 1, 


Unit B3, Theorem 3.2 
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and hence obtain 


Thus J ^ is of the form given in Condition (3.5). 


i _ Jl 

z z3 


for |z| > 2. 


It follows, by the Flow Mapping Theorem with a = 1, that the velocity function 

qiz)=qi4J-\z))iJ-^)'{z) 


solves the Obstacle Problem for K with circulation 27rc ziroimd K. To simplify 
this formula, we note that, by the Inverse Function Rule, 

^ l-l/(J-i(z))2’ 

and hence 


J{w) = w + 1/w 
J'{w) = 1 - 1/w^. 


\j-^iz)-l/J-^iz)J 


—(vAtp 

as required. ■ 

Remarks 

1 Notice that if c = 0 in Example 3.1, then we obtain the expected solution 
q{z) = 1, given in Problem 3.3. 

2 Note that we can compute the velocities on either side of the obstacle 

K = [-2,2], by putting z = x + ly in the expression for q{z) and letting y tend 
to 0 through positive values (written y —> O^) and through negative values 
(written y —» 0“). 


For example, put x = 0, so that z = iy. Then 

z / i - 4/z2 = / fory>0, 

\ -j|y|\/l -4/(-y2), fory<0, 

{ Wy^ + 4, for y > 0, 

+4, for y < 0, 

r 2i as y —* O’*", 

^ \ -2i as y —*0~. 


{ 1 — ^c as y -+ 0+, 
1 + ^c as y -+ 0~. 


This suggests that negative circulation reinforces the flow above the obstacle 
and opposes it below. 


w = J-\z) 

z = J{w) 

- w+ - 

= J-\z) + \/J-\z). 
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3 An alternative method of obtaining the velocity function in Example 3.1 is 
first to determine the complex potential function f 2 = ^la,c ° f (given by 
Theorem 3.2), differentiate f2 and then take its conjugate to give the rule for 
q — Q'. We ask you to carry out this approach in the next problem. 


Problem 3.5___ 

Solve the Obstacle Problem for K = [-2,2], with circulation 27rc around K, by 
first showing that a complex potential function for the solution is 

0(0) = z — zcLog^|(z + z-\/l — 4 /z^) \ . Note that if c = 0 then 

n(.z) = 2. This happens because 

{Hint: The identity = J - fcLog will help.) , n 1 

_ Sli.ofz) = 2 + - = J( 2 >) 

and 

The complex potential function fl found in Problem 3.5 should, in principle, /(^) = 

enable us to plot the streamlines for the corresponding flow round K = [—2,2], 
since these have equations of the form 

Im(f2(2)) = constant. 


However, if c 7 ^ 0, then this equation is so complicated that the streamhnes can 
only be plotted accurately with the help of a computer. 


An alternative approach to sketching these streamlines is to use the fact that 
the conformal mappings f and f~^ map streamlines to streamlines (see 
Remark 1 after the Flow Mapping Theorem). Thus, if P is any streamline for 
the velocity function qi^c, then f~^{T) is a streamline for the velocity function 
q. We can therefore use the function f~^ = J to transfer the appropriate 
streamline diagram in Figure 3.2 (where we had c < 0) from C - ATi toC — K. 
Of course, you can only hope to do this approximately by hand, but some 
precision can be achieved if you first calculate the effect that f~^ = J has on 
the stagnation points of the flow with velocity function qi^c (which were found 
in Problem 3.2). Since a = 1 in this case, the stagnation points are 

w= lic± 0 ^^^, if -2 < c < 0, (3.7) 

and the single point 



To illustrate the procedure we consider the case -2 < c < 0, for which the 
stagnation points lie on the boundary dKi = {w : jiyj = 1} of K^. The function 
= J IS one-one and conformal on the region C - iPi, but it can be 
extended to dKi by using the rule J{w) = w-\-l/w there. 


The material from here up to 
and including Figure 3.8 may be 
omitted on a first reading. 


It can be shown that J maps both the upper half and the lower half of dK\ 
onto K = [—2,2], and so the restriction of J to dKi is not one-one. Also, it can 
be shown that J maps points in C - ATi near the upper half of dKi to points 
just above K and points in C — iPi near the lower half of dKi to points just 
below K (see Figure 3.7). 



Figure 3.7 


Since the stagnation points given by Equation (3.7) lie on the lower half of 
dKi , the streamlines of the flow which end at these two stationary points are 
mapped by /“^ = J to streamlines of the flow with velocity function q which 
approach K = [—2,2] from below, ending at the points j(^^ic± j. 
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Once these points have been determined, we can sketch the corresponding 
streamlines (and others) by using the fact that 

q{z) = 1, for \z\ large, 

together with the symmetry of the streamhnes with respect to the imaginary 
axis, which follows from the corresponding symmetry of the streamlines in 
Figure 3.2(b) and that of the obstacle K. A resulting sketch (for c = -\/3) is 
shown in Figure 3.8. (The streamlines corresponding to other obstacles K may 
be sketched in a similar way.) 



Figure 3.8 


To increase the range of obstacles which we can deal with, we now introduce 
the family of mapping functions 

Ja{z)^z+^ (2 eC-{0 }), 

where a > 0. These functions are closely related to the Joukowski function 
J = Ji, and we can deduce their basic properties by observing that 

J.lz) = <i(z/a+^y. 

that is, w = Ja{z) is obtained by composing 

(a) the scaling Zi — zja-, 

(b) the Joukowski function Z2 = J{z\) = zi + 1/zi; 

(c) the scaling w = az2- 

Problem 3.6- 

(a) Show that Ja is a one-one conformal mapping from C- Ka onto 
C - [-2a, 2a], and that Ja maps dKa onto [-2a, 2a]. 

(b) Verify that the inverse function of the restriction of Jo to C - Ka is 

= \{w + (w e C - [-2a,2a]). 


The results of Problem 3.6 show that, if we reverse the roles of z and to, then 
we can use J“^ as the function / in the Flow Mapping Theorem to solve the 
Obstacle Problem for K = [-2a, 2a], where a > 0. 

Problem 3.7- 

Use the function J~^ to show that the solution to the Obstacle Problem for 
K = [-2a, 2a] with circulation 27rc around K is 


Note that 

Ua.c = Ja-ic Log. 


Note that Problem 3.7 
generalizes Example 3.1. 


The functions Ja can also be used to solve the Obstacle Problem for flow 
around various other obstacles, such as an ellipse. 






Problem 3.8___ 

(a) Show that if r > a, then Jq maps the circle Cr = {w : |u;| = r} onto the 
ellipse in the z-plane shown in Figure 3.9. 





(b) Deduce from part (a) that the solution to the Obstacle Problem for the 
obstacle K comprising the elhpse in part (a) together with its inside, with 
circulation 27rc around K, is 


q{z) = 1 - 


( r^-a^ + icJa\z) \ 


where J~^(z) = ^(z + Zy/1 - 4a^/z^y 

{Hint: Use the solution to Problem 3.7 when estabhshing 
Condition (3.5).) 


We end this subsection with a proof of the Flow Mapping Theorem. This proof may be omitted on a 

first reading. 

Theorem 3.2 Flow Mapping Theorem 

Let AT be an obstacle and let / be a one-one conformal mapping from 
C - K onto C — Ka, where a > 0, such that 

f{z) =z + ao + ^ + ^ + ---, for |z| > A, (3.5) 

where R> 0 and oq, a_i, a_2,... 6 C. Then the velocity function 

9(z) = 9 o,c(/(z))7^ (zgC-K) (3.6) 

is the unique solution to the Obstacle Problem for K, with complex 
potential function 

D = fia,c ° /■ 


Proof Before verifying Properties (a), (b) and (c) of the Obstacle Problem, We omit the proof of the 
we observe that uniqueness of the solution q. 

9(^) = 9o,c(/(2))7^ 


is a model flow velocity function on 71 = C — K, since it is the conjugate of an 
analytic function there. 


Now, using the Laurent series (3.5), we have 

Um = 1 and hm f'{z) = 1. 


Note that the functions 
z f{z)/z 

and /' both have removable 
singularities at oo; see Unit Dl, 
Section 1. 
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It follows that 


t 

I 


\ 


e 


lim q{z) — lim ( 1 - 


if{z)r 


f{^)J 


m 


which verifies Property (a) of the Obstacle Problem. 

To verify Property (b) of the Obstacle Problem, we use the fact that / is 
one-one and analytic on V, so that, by the Inverse Function Rule, / has an 
inverse function f~^ which is analytic on C — Since is analytic on Cm 
it follows that the composite function f2a,c ° / is analytic on the region 

r'(C.-R'a) = 7^-r'(Sa). 

where Sq = {u € R : u < —a} is that part of the negative real axis in C — /fa 
(see Figure 3.10). 




Since Eq is a simple smooth path and / ^ is one-one and analytic, it follows 
that E = f~^{T,a) is also a simple smooth path. Moreover Q = fla,c o/ is a 
complex potential function for 9 on — E, because 

(Chain Rule) 


= q{z), for 2 e 7^ — E. 

To complete the verification of Property (b), we need to show that there is a 
constant k such that 

for each a £ dK, we have lim Im(fl(2)) = k. (3.8) 

To do this, first note that if a € dK, then a is a limit point of 7^ - E; indeed, 
each open disc with centre a must contain a point of 7^ — E (because a € dTZ 
and E cannot fill up an open set) and so we can construct a sequence {z„} in 
7^ - E such that Zn a. 


Now suppose that {zn} lies in 7^ - E and —> a. Then the sequence 

{•Wn} = {f{zn)} lies in C,r — Ka, and we claim that |u;„| —» a. Indeed, if e > 0 
is given ajid Fj = {ly : |w| = a -I- e}, then is a simple-closed contour in 

TZ which surrounds, but does not meet K (see Figure 3.11). Moreover, f~^ 
maps the outside of Fe onto the outside of /"^(F^). It follows that the terms of 
{zn} must eventually lie inside /“^(Fj), and hence that the terms of {wn} 
must eventually lie between F^ and dKa- Hence —» a as n —» 00. 



Figure 3.11 


These plausible 2issertions can 
be verified by using the 
Argument Principle (Unit C2, 
Theorem 2.2), together with the 
fact that / is a one-one analytic 
function of the form given by 
Equation (3.5). We omit the 
details. 
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Therefore, if u;„ = u„ + iVn, for n = 1,2,..., then we have 
Im(n(^„)) = Im(n„,c(/(^„))) 

= Im(no,c(iyn)) 

= Un (l - - clogg |u;„| (by Equation (3.3)) 

—cloggO as n —» oo. 

Thus Condition (3.8) holds with k = —clog^a. 

Finally, we check Property (c). If T is any simple-closed contour 
surrounding K, then /(T) is a simple-closed contour surrounding Ka and so 

Cr + = J dz 

= J^^{f{z))f'iz)dz 

= / q^{w)dw = 2nc (by Theorem 3.1(c)). 

Jf(r) 

Hence Cr = 27rc, as required. ■ 


The substitution 

^ = /(«), dw = f'{z)dz 
can be justified as in the proof 
of the Argument Principle. 

Note that this integral also gives 
J"r = 0. 


3.4 Aerofoils 

You saw in Problem 3.6(a) that the analytic function Ja{w) — 'W + a'^/w maps 
the circle Ca = {w : |u;| = a} to the line segment [—2a,2a], as well as being a 
one-one mapping of the region exterior to Ca onto C — [—2a, 2a]. We consider 
next why the image under Ja of the smooth path Ca is ‘pointed’ (non-smooth) 
at both ends. This leads on to the definition of a type of aerofoil, which is a 
shape pointed at just one end. 

As shown in Unit A4, if / is an analytic function aud a is a point of its domain Unit A4, Theorem 4.2 

at which /'(a) 0, then / is conformal at a, which means that / preserves the 

angle between two smooth paths which meet at a. However, if a is a point for 

which f'{a) = 0 but /"(a) ^ 0, then the effect of / is to double the angle 

between two smooth paths emerging from a. (For a justification of this 

statement, see the Local Mapping Theorem {Unit C2, Theorem 3.2), and the 

Remark following its proof in Subsection 3.3, taking n = 2.) We call a point a 

at which the analytic function / has zero derivative a critical point of /. 

This explains why Ja(C'a) = [-2a, 2a] is ‘pointed’ at both ends. For 

J'aiw) = 1 - ^ 

is zero at ta = ±a, but J" (±a) ^ 0. The image points of lu = ±a under Ja are 
at z = ±2a. The circle Ca may be thought of as a pair of smooth paths 
meeting at the points w = ia, with an angle n between them. The images of 
these two paths therefore meet at an angle 27r, producing the pointedness of 
the image contour. The image of any smooth path passing through w = ±a has 
a similar pointed nature. 
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Problems 3.6(a) and 3.8(a) 


Figure 3.12 illustrates two of the results obtained so fax. The function 
Jaiw) =w + a^/w maps 

(a) the circle Ca = {w: |ii;| = a} to the line segment [—2a,2a]; 

(b) the circle Cr = {w: |r(;| = r}, where r > a, to the ellipse 



Figure 3.12 


Consider now the circle 

C = {w ; |w + 6| = a + 6}, 

where b = |(t- - a). This circle, which is also shown in Figure 3.12, has the 
same tangent as the circle Ca at w = a, so that its image JaiC) can be 
expected to resemble that of Ja{Ca) = [—2a, 2a] in the vicinity of 0 = 2a. In 
particular, Ja{C) should be ‘pointed’ at z = 2a. The circle C has the same 
tangent as the circle Cr at w = —r, so that its image JaiC) can be expected to 
resemble that of the ellipse JaiCr) in the vicinity oi z = —r — a?/r. 

At all points of C other than w = a we have J'aiw) ^ 0, so the image contour is 
smooth except at the single point z = 2a. Also, the mapping Ja preserves the 
symmetry of C about the real axis, since Jaiw) = Jaiw). 

These considerations lead to the sketch of Ja{C) shown in the z-plane on the 
left of Figure 3.12, and repeated in Figure 3.13. 



Figure 3.13 


The obstacle with boundary Ja{C) is called a Joukowski aerofoil. Note that the 
circle C which defines it passes through the critical point w = a of Ja, whereas 
it surrounds the other critical point w = —a. 


In order for a Joukowski aerofoil 
to be symmetric, the centre of 
the circle C must lie on the 
(negative) real axis, as in 
Figure 3.13. Otherwise, the 
resulting aerofoil is asymmetric. 


Definition A Joukowski aerofoil is an obstacle which (possibly after 
an appropriate translation or rotation) has boundary JaiB), where B is a 
circle which passes through the critical point w = a of Ja, and surrounds 
the other critical point w = —a. 

The point z = 2a is called the trailing edge of the aerofoil. 


Remark The boundary JaiB) is smooth except at z = 2a, and surrounds the 
point z = —2a. 





Problem 3.9_ 

(a) Suppose that C is the circle {w ; |w + 6| = a + 6} (where a,b> 0), and that 
K is the Joukowski aerofoil for which dK = Ja{C). Show that the solution 
to the Obstacle Problem for K with circulation 27rc around K corresponds 
to the complex potential function 

n(i:) = w + b+ -icLog(u; + 6), 

where w = J~^{z). 

(b) Determine, in terms of w, the flow velocity function q for this flow, and 
show that ^liin^ q{z) (the limiting velocity at the traihng edge) exists if and 
only if c = 0. 


Apart from a case involving translation (Problem 3.4), we have so far 
considered only examples of obstacles which are symmetric about the real axis. 
While an aeroplane is in flight, the wing cross-section is at a shght angle to the 
oncoming airstream. We may model this situation by placing a symmetric 
Joukowski aerofoil at an angle 0 to a uniform stream, and attempt to And a 
complex potential for this version of the Obstacle Problem (see Figure 3.14). In 
these circumstances, 0 is known as the aerofoil’s angle of attack. 



In order to make progress here, we shall extend the family of mapping 
functions Jq. The extension required is to allow the parameter a to become 
complex, that is, to consider the family of mapping functions 


Ja{w) =w - 


(3.9) 


where a is any non-zero complex number. Note that if 0 = - Arg a, then 
Equation (3.9) can be written as 

showing that the transformation z = Ja{w) is a composition of the three 
transformations 

(a) wi = (an anticlockwise rotation through the angle 0); 

(b) W2 = 7|a|(u;i) = lui -h lap/ioi (a mapping of the type Jo, with a = |a|); 

(c) z = W2e~^’^ (a clockwise rotation through the angle 0). 


Thus, if we denote by R^j, the anticlockwise rotation of C about 0 through the 
angle 0, then 

Jc = o J|„| o R^. (3.10) 

Suppose that K\a\ is the closed disc {w: |i/;| < |a|}, whose boundary 
C'laj = includes the point a. Figure 3.15 shows, using the decomposition 

of Equation (3.10), that Ja maps C — A'|q,| onto C — L{—2a, 2a), as well as 
mapping C\a\ to L(—2a,2a). 


In other words, we have 



where —tt <<!><%. 


L{—2a, 2a) denotes the line 
segment joining —2a and 2a. 
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Figure 3.15 

Since each of the constituent mappings is one-one on the appropriate region, 
the mapping Ja restricted to C - K\a\ is also one-one. This restriction of Ja 
therefore has an inverse function which is a one-one conformal mapping 
from C — L(-2Q,2a) onto C - From Equation (3.10), we have 



Using the formula for ^{z) (see Problem 3.6(b)), with a = |q|, we deduce 
that 



that is, 

= + j (zeC-X(-2a,2a)), (3.11) 

as expected in view of the formula for J~^. 

Example 3.2 

(a) Find a complex potential which provides a solution to the Obstacle 
Problem when the obstacle K is the line segment L{—2a,2a) and the 
circulation around K is 27rc. 

(b) Verify, from your answer to part (a) with c = 0, that the circulation-free 
complex potential for a = a (where a is reed) is Q{z) = z, as expected. 

(c) Determine the velocity function q which solves the Obstacle Problem for 
K = L{—ai,ai), with zero circulation. 

Solution 

(a) As shown above, J~^ is a one-one conformal mapping from C — L(—2a, 2a) 
onto C - K\a\- Also, it can be shown that the function / = J~^ satisfies 
Condition (3.5) of the Flow Mapping Theorem. The appropriate Laurent 
series is obtained as in Problem 3.7, with a in place of a. A suitable 
complex potential for the flow past K, with circulation 27rc around K, is 

J^(^)=(n|a|,cOJ-')(0) 

= JaHz) + - icLog{j-\z)) 

= ^ - icLog(J-^(z)), 



Here K may be regarded as a 
thin plate of height 2a placed at 
right angles to a uniform flow. 
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(b) With c = 0, we have 



so that if a = a, where a is real, then the complex potential is 

n(z) = z, 

as in the solution to Problem 3.3. 

(c) If c = 0 and a = ^ai, then we find from part (a) that 




The required flow velocity function is therefore 


q{z) = n'{z) 




— 1 /. ■ Note that q is unbounded near 

±ai, the ends of the line 

Remark in part (c), the velocity function q could have been obtained directly segment, 
from Equation (3.6), as we did in the solution to Example 3.1. 

To end the section, we ask you to solve the Obstacle Problem for a rotated aerofoil. 


Problem 3.10_ 

(a) Suppose that the Joukowski aerofoil K shown in Figure 3.14 is obtained by 
rotating the aerofoil in Figure 3.13 and Problem 3.9 clockwise about the 
origin through an angle (j). We then have 

dK = (R-^oJ^) (C), 

where, as before, C is the circle {lu : |tt) + 6| = a + 6}. With a = use 

Equation (3.10) to show that the aerofoil boundary dK is the image under 
Ja of the circle 


R^\C)={w:\w + be-'*\=a + b}, 



Figure 3.16 
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(b) Determine a complex potential which provides a solution to the Obstacle 
Problem when the obstacle K is the Joukowski aerofoil of part (a), at angle 
<j) to the incoming stream, with circulation 27rc around K. 

(c) Determine the corresponding flow velocity function q, in terms of 

w = and show that lim q{z) (the limiting velocity at the traihng 

edge) exists if and only if 
c = -2(a + b)sm<(). 

{Hint: In order to make the algebra more manageable, use the auxiliary 
variables 

w\=w — a = w — ae~*^, 

1^2 = 10 + be~''^ = loi + (a + b)e~'^, 
and note that ici -+ 0 as z —» 2a.) 


The consideration of aerofoils is followed up in the next section, where we 
develop our model further in order to calculate the force acting on an obstacle, 
due to the fluid streaming past it. You will see that the circulation around the 
obstacle is of primary importance in this context and, as far as aerofoils are 
concerned, the solution to Problem 3.10(c) will play a key role. 



4 LIFT AND DRAG 


After working through this section, you should be able to: 

(a) explain the meaning of Bernoulli’s Equation, and apply it to find the 
pressure distribution of a model fluid flow; 

(b) find the total force acting on an obstacle due to a model fluid flow, either 
directly from the pressure distribution, or by the application of Blasius’ 
Theorem, or by using the Kutta-Joukowski Lift Theorem; 

(c) understand why the Kutta-Joukowski Hypothesis leads to a bounded 
velocity function for flow past a Joukowski aerofoil, and a unique value for 
the force on the aerofoil. 


4.1 Bernoulli’s Equation 


Our aim in this section is to investigate the force acting on an obstacle located 
in a fluid flow, using the velocity functions found in various cases in Section 3. 
This force is caused by the fluid pressure on the boundary of the obstacle, and 
its value may be obtained by calculating a suitable contour integral. 

Pressure is a scalar quantity (that is, it has magnitude but no associated 
direction), which may be modelled within the flow region by a non-negative 
real-valued function p. This function is called the pressure distribution within 
the fluid. It is defined to be the magnitude of the force per unit area acting 
within the fluid. This means that, on any small area 6A of a fixed surface 
which includes z, the pressure gives rise to a force of approximate magnitude 
p{z)6A, in the direction perpendicular to the surface at 2 (see Figure 4.1). 

We now relate the pressure p{z) of the fluid at the point 2 to its velocity q{z). 
A convenient modelling link is provided by the equation below. 


Force is a vector quantity. 

The pressure on the surface of a 
solid object in contact with a 
fluid is associated with 
momentum changes brought 
about by large numbers of fluid 
molecules ‘bouncing off’ the 
molecules in the solid surface. 
This explanation lies outside our 
continuum model for the fluid. 


v{z)6A 


Figure 4.1 


Bernoulli’s Equation 

In a fluid of constant density p, the relationship between the fluid 
pressure p{z) and the fluid speed |g(2)|, at the point 2, is given by 
Bernoulli’s Equation 

P{z) + \p\<l{^)? = Po, (4.1) 

where po is a constant. 


Remarks 

1 A derivation of Bernoulli’s Equation is given in Subsection 4.4. It is an 
integral version of Newton’s Second Law (force = mass x acceleration) applied 
to a small volume of fluid, and is the counterpart of an energy equation in 
single-particle mechanics. Its derivation assumes that all forces acting within 
the fluid other than those due to the pressure may be ignored. In particular, 
gravitational and viscous effects are ignored. 

2 The constant po is the fluid pressure at any stagnation point, as can be seen 
by putting q{z) = 0 in Bernoulli’s Equation. For a flow satisfying the Obstacle 
Problem, po can be expressed in terms of poo, the pressure at large distances 
from the obstacle. Since lim q{z) = 1, we have 


Daniel Bernoulli (1700-82) came 
from a Swiss family which 
yielded famous mathematicians 
throughout the eighteenth 
century. He was in succession 
Professor of Medicine, 
Metaphysics and Natural 
Philosophy at the University of 
Basle. The word 
‘hydrodynamics’ originated from 
his publication of 
Hydrodynamica in 1738, which 
was the first substantial work in 
the theory of fluid motion. 


Viscous effects are caused by the 
internal friction of a fluid. We 
say more about this in 
Subsection 4.3. 


Po = Poo + \p- 

3 In words, Bernoulli’s Equation says that a pressure increase is accompanied 
by a decrease in the fluid speed, and vice versa. 
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The prediction in Remark 3 may be used in circumstances which go beyond our 
current two-dimensional model. For example, if you hold two sheets of paper 
vertically, a few centimetres apart, and blow between them, you will find that 
the sheets tend to move together (see Figure 4.2). The explanation is as follows. 

The blowing between the paper sheets produces a higher flow velocity in the 
region between them than outside. Bernoulli’s Equation predicts that the 
pressure between the sheets will therefore be less than that outside, and this 
pressure difference provides a force on eeich sheet, pushing it towards the 
centre. The sheets therefore move together. 

Another example of Bernoulh’s Equation in action requires a little more by way 
of equipment, namely, a table-tennis ball and a thin tube. You will find that 
the table-tennis ball can be balanced on the vertical airstream created by 
blowing through the tube, as illustrated in Figure 4.3. The ball’s situation is 
stable (for as long as you can keep blowing at a steady rate). If pushed a small 
distance sideways and then released, it will return to the centre of the stream. 



\I/ 


Figure ^,3 I blow 


Problem 4.1_ 

Use BernouUi’s Equation to explain the stability of the table-tennis ball just 
described. 


Having invited you to convince yourself that the qualitative predictions of 
Bernoulli’s Equation are reasonable, we return to our two-dimensional flow 
model, and in particular to the flow past a disc of radius a, considered in 
Subsection 3.1. The velocity function in this case (with circulation 27rc around 
the disc) is 

9o,c(^) ~ ^ ~ ^ ~ 


Problem 4.2_ 

(a) Use Bernoulli’s Equation to show that the pressure distribution of this flow 
past the disc is 


P{z) =Po- 



(b) Find the pressure distribution on the boundary of the disc, 
{ae*‘: 0 < t < 27r} , 


giving your answer in terms of po, p, c/a and sint. 

(c) Deduce that the pressure distribution on the boundary is symmetric about 
the imaginary axis, and that if c = 0 it is also symmetric about the real 
axis. (As a result, the horizontal component of the total force on the 
obstacle due to the fluid is zero, and the vertical component is also zero if 
c = 0.) 


sheets held 





Figure 4-^ 


A similar effect may be observed 
by placing a (not too ripe!) 
tomato, or other round object, 
in a kitchen sink beneath a 
steady stream of water from the 


Equation (3.4) 
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Prom an expression for the pressure distribution on the boundary of an 
obstacle, as obtained in one particular case in Problem 4.2(b), it is possible to 
deduce the total force caused by pressure exerted on the obstacle due to the 
fluid flow. We show how to obtain this force in the next subsection. 


4.2 Total force on the obstacle 


The results in this subsection establish the total force on any obstacle caused 
by a flow which is a solution to the Obstacle Problem. This total force is 
caused by the effects of pressure around the boundary of the obstacle. Pressure 
is magnitude of force per unit area, but in our two-dimensional model one of 
the dimensions of area is directed into the page, perpendicular to the complex 
plane which we are using. Hence, the effect of the pressure distribution p, on a 
small segment of the (cross-sectional) boundary of length 6s, is a force (per 
unit length into the page) of approximate magnitude p(z)6s, where z lies on 
the boundary segment concerned. This force is directed into the obstacle and 
perpendicular to the boundary at z, as illustrated in Figure 4.4. 



The total force on the obstacle K is the ‘sum’ of all such contributions around 
the boundary dK, where the sum has to take account of the directions as well 
as the magnitudes of the force contributions. 

In order to calculate this force we now suppose that the boundary dK of the 
obstacle is a simple-closed smooth path of length L, which is traversed 
anticlockwise with unit-speed parametrization 7: [0, L] —> dK (so that 
7(L) = 7(0)). Then, for each s € [0,1-], the unit vector 7'(s) is a tangent vector 
to the path and directed in the anticlockwise sense, as shown in Figure 4.5. 

This figure also shows that the effect on 7'(s) of multiplying by i = is to 
rotate the vector by 7r/2 anticlockwise, producing another unit vector i'y'{s), 
which is directed into K and perpendicular to the boundary dK at the point 
z = 7(s). The corresponding force contribution (directed along i'y'{s)), for a 
length 6s of the boundary which includes 7(5), is approximately 



7'(«) 

Figure 4.5 


p(7(s))<5s X i-y'{s) = ip{-y{s))'y'{s)6s. 


From these observations, it can be shown that the total force F acting on K 
(per unit length into the page) is 


( 4 - 2 ) 

Special names are given to the real and imaginary parts of the force F. 


Definitions if F is the total force acting on an obstacle, due to a fluid 
flow whose velocity function satisfles the Obstacle Problem, then 

(a) Re F is called the drag on the obstacle; 

(b) Im F is called the lift on the obstacle. 


Note that F is a (complex) 
constant. With an appropriate 
change to the limits of 
integration, the same formula is 
valid for any equivalent 
parametrization of dK. 


For the boundary of the disc of radius a, considered in Problem 4.2, we have 
the unit-speed parametrization 

7(s) = (s e [0,27ra]), 

which is obtained from the standard parametrization 7i(f) = ae*‘ {t G [0,27r]) 
by putting s = at. Now 

7'(s) = = ze’*. 
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and so, in this case, Equation (4.2) becomes 


/•2’f 

F = i j p[ae'^) ie'^adt 
= —a j p(ae'^) e'^dt. 


Problem 4.3- 

By using the above formula for F, aud the expression 
p(ae‘‘) =po - 1(2 Sint-, 

from the solution to Problem 4.2(b), find the total force (per unit length into 
the page) on a cylinder of radius a, due to a uniform flow with circulation 27rc 
around the obstacle. What are the lift and the drag in this case? 


The result which you found in the particular case in Problem 4.3 (that is, zero 
drag and positive lift for negative circulation) generalizes to edl obstacles. In 
order to show this we first re-express Equation (4.2) as follows. 


Theorem 4.1 Blasius’ Theorem 

Let q be a solution to the Obstacle Problem for an obstacle K, whose 
bound2iry dK is a simple-closed smooth path, and suppose that q has a 
continuous extension to dK. Then the conjugate of the total force acting 
on the obstacle is 

F = ^ip [ {q{z))^ dz, (4.3) 

JdK 

where p is the density of the fluid. 


This result was published by 
H. Blasius in 1910. A similar 
theorem, which also bears his 
name, states that the moment 
(turning effect) M about 0, 
experienced by the obstacle due 
to the flow, is given by 



Remark if K is an obstacle whose boundary is not a simple-closed smooth 
path, then the integral in Equation (4.3) may not be defined. In this case we 
take the conjugate of the force acting on the obstacle to be 

F = ^ip (4.4) 

where P is any simple-closed contour surrounding K. This is plausible for two 
reasons; 

1. the integral in Equation (4.4) is independent of the choice of simple-closed 
contour P surrounding K because the function z i —* (9(2))^ is analytic on 
C — K, and so Cauchy’s Theorem can be apphed as in the proof of the 
Shrinking Contom: Theorem {Unit B2, Theorem 1.4); 

2. it can be shown that if dK is a simple-closed smooth path, then the 
integrals in Equations (4.3) and (4.4) are equal. 


Proof On using Bernoulh’s Equation to substitute for p(7(s)) in 
Equation (4.2), we obtain 

This may be simplified by noting that 

f Po7'(s) ds = poi-yiL) - 7(0)) = 0. 

Jo 
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The complex conjugate of F is therefore 


= k( 7 (s))P 7 '(s)rfsj 

= k(7(s))p7'(s)<is. 

Now since 7 is a unit-speed paxametrization, we have 7'(s) = e*"^, for some (t>, 
and hence 7'(s) = 

Also, dK is made up of streamhnes of the flow, so that if qris) denotes the 
component of 9(7(5)) in the direction of 7'(s), then 

9(7(5)) = 9T(s)e*'^ and 9(7(5)) = 97(5)6“*'^. 

Hence 

19(7(5))!^ = 9 ( 7 (s)) 9 ( 7 (s)) = ( 9 t( 5 ))^ = ( 9 ( 7 ( 5 )))=^e 2 ‘<^, 
and so the integrand is 

|9(7(5))P7'(5) = = (q(7(s)))^j'(s). 

Thus we obtain the required formula 

[ (9(7(s))f7'(s)ds = \ip f {q{z)fdz. ■ 

Jo JdK 

Problem 4.4_ 

Show that applying Blasius’ Theorem to the conjugate velocity function 

_ , ic 

9a,c (^) = 1 ^ - 7 

gives a more rapid derivation of the result of Problem 4 . 3 . 


The next theorem is the promised generalization of the result of Problem 4 . 3 . 

Theorem 4.2 Kutta-Joukowski Lift Theorem 

Let 9 be a solution to the Obstacle Problem for an obstacle K, with 
circulation 27 rc aroimd K. Then the obstacle experiences a total force 

F = — 2 ‘Kcpi, 
that is, 

(a) a zero drag, and 

(b) a lift of — 2 ncp. 


Proof Since the velocity function 9 is a solution of the Obstacle Problem, its 
conjugate 9 is analytic on C - A and lim q{z) = lim q{z) = 1. Hence 9 has a 
removable singularity at 00 and so, if K is contained in the disc {z : \z\ < R}, 
then 9 has a Laurent series of the form 

q{z) = 14- -I- -I-, for |z| > R. ( 4 . 5 ) 

We can now express a_i in terms of the circulation 27 rc around K. Prom 
Equation ( 1 . 7 ), we have 

Cr + iFt = j q{z) dz, 

where Cr and Ft: 2ire, respectively, the circulation along and flux across emy 
simple-closed contour T which surrounds K. 


Wilhelm Martin Kutta 
(1867-1944) spent most of his 
life in Munich, where he became 
a professor. He was a pioneer of 
two-dimensional wing theory, 
and discovered this result in 
1902. It was found 
independently by Joukowski in 
1906. 


This proof is for the case where 
dK is a simple-closed contour to 
which 9 has a continuous 
extension. 

The function g{w) = q{l/w) has 
a removable singularity at 0, 
since 

lim g{w) = 1. 
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Talcing F to lie outside the disc {z : \z\ < R}, we obtain 

= 2ma-i. 

The integral J {^ + H-^ dz is zero because the integrand has a 

primitive on {z : |z| > R}; see Step (e) of the proof of Laurent’s Theorem in 
Unit B4, for similar reasoning. 

Since Jr = 0 (see Remark 4 following the Obstacle Problem) and Cp = 27rc, we 
deduce that 


a_i = -ic. (4.6) 

Now, by the remark following Blasius’ Theorem, 

F= \ip j^{q{z))^dz, 

where, as above, the simple-closed contour F lies outside the disc {z : \z\ < R}. 
Hence, using Equations (4.5) and (4.6), we have 

= ^ip{—2ic)2m = 2-Kcpi, 
so that F = —2ncpi, as required. ■ 

With this result in mind, we return to the subject of aerofoils. In 
Problem 3.9(b), you saw that a symmetric Joukowski aerofoil with its axis on 
the real axis can have a finite limiting velocity lim q{z) at its traihng edge if 
and only if there is no circulation. In Problem 3.10(c), this was generalized to 
the case of a symmetric aerofoil with angle of attack 4), showing that hm q{z) 
exists if and only if 

c = —2(a -I- b)sm<t>, 

where 27rc is the circulation around the aerofoil. 

In each case, if the circulation does not have the specific value indicated, then 
the flow velocity q is imbounded near the traihng edge. Such imboimded values 
in the model do not correspond with the real world, and this provides a 
motivation for avoiding such a situation. Moreover, Bernoulli’s Equation 
predicts that negative pressure will be encountered in a region where q is 
unbounded. However, negative pressure is not an effect observed in nature. 

The only way around this difficulty, within the scope of our current model, is to 
assume that the circulation takes precisely the value required to avoid the 
problem. This is the purpose of the following hypothesis. 


Kutta-Joukowski Hypothesis 

The circulation around a Joukowski aerofoil is such that the flow velocity 
q is bounded throughout the fiow region, and q{z) tends to a hmit as z 
approaches the traihng edge of the aerofoil. 
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Remarks 

1 This hypothesis, as well as being convenient in pinning down a specific veilue 
for the circulation, is in good agreement with experimental observations. 

2 For a value of c which satisfies the Kutta-Joukowski Hypothesis, there is 
just one stagnation point on the boundary dK, whereas for nearby values of c 
there are two stagnation points. Roughly speaking, therefore, the effect 
described by the hypothesis is achieved by ‘placing a stagnation point (zero) on 
top of a singularity’. 

Problem 4.5_ 

How does the lift on a symmetric Joukowski aerofoil depend on the angle of 
attack 4), assmning that the Kutta-Joukowski Hypothesis holds? Does this 
agree with what you know of aeroplane wings? 


The solution to Problem 4.5 suggests that the model which we have 
constructed is only partially successful in describing fluid flows in the real 
world. In the next subsection, we consider just how good the model is. 


4.3 Successes and limitations of the model 

In this subsection, we assess the usefulness of studying model fluid flows. Prom 
the point of view of the current course, this unit has demonstrated the 
applicability of many of the complex analysis results and ideas introduced 
earlier. While this should have helped you to understand and appreciate these 
ideas further, it is not a reason to accept the model imcritically. 

On the plus side, the streamline pictures for the solution to the disc version of 
the Obstacle Problem, for various values of the circulation (see Figure 3.2), can 
be shown experimentally to have some validity, and a lift due to the circulation 
is observed. Also, for small values of the angle of attack presented by an 
aerofoil to the oncoming airstream, the prediction of lift being increased by 
raising the leading edge of the wing (Problem 4.5) is certainly borne out. In 
level flight, this upward lift force balances the downward pull due to gravity. 
Vaxiation of the angle of attack allows an aircraft to chmb (greater angle than 
for level flight) or descend (lesser angle) in a controlled fashion. 

There is, however, a definite hmit to the validity of this prediction, as pointed 
out in the solution to Problem 4.5. If the angle of attack becomes too laxge 
then staffing occurs; that is, the strength of the lift force decreases dramatically, 
and an aeroplane connected to the aerofoil would fall out of the sky. 

There is also the question of the physical mechanism which produces the 
circulation necessary for lift to be provided. Our model contains circulation 
almost as an arbitrary ingredient whose value can be altered at will, except in 
the case of an aerofoil where the Kutta-Joukowski Hypothesis provides some 
rationale for picking a specific value for circulation. Even this rationale may 
seem dubious, since it arises from the desire to eliminate unbounded velocity 
values, though there is good experimental support for it. 

The Kutta-Joukowski Lift Theorem is curious in two respects. Firstly, it 
predicts zero drag, that is, an absence of any force component in the direction 
of the incoming fluid stream. This prediction is seriously at odds with 
observation: for example, the table-tennis ball in a jet of air (see Problem 4.1 
and the text before it) is maintained by the jet at some distance above the top 
of the tube through which air is blown, at the point where the drag due to the 
jet exactly balances the weight of (gravitational force on) the ball. In the 
absence of drag, the ball would descend to the mouth of the tube. 

Secondly, the Kutta-Joukowski Lift Theorem predicts, for any obstacle, a lift 
force which depends only on the circulation and fluid density, and is therefore 


This subsection is intended for 
reading only. 


The mismatch here between the 
prediction of the model and 
experimental observation is 
known as d’Alembert’s Paradox. 
Jean Le Rond d’Alembert 
(1717-83) was a leading French 
mathematician, a philosopher 
and a friend of Voltaire. His 
Traite de I’Equilibre et du 
Mouvement des Fluides was 
published in 1744, and the 
Theorie Generale des Vents in 
1745. 
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independent of the object’s shape or dimensions. This seems to suggest that 
the aerofoil shape is unnecessary, and that any obstacle, no matter what its 
profile, would be lifted equally by a given circulation. Again, this prediction is 
not what is observed in practice, although the circulation can depend to some 
extent on the shape of the body itself. 

Clearly our model is lacking in some significant respect. Its major 
shortcomings arise because we have assumed the fluid flow to be locally 
circulation-free. In fact, non-zero circulation around small closed contours is 
always present in a moving fluid sufficiently close to a solid boundary. The 
mechanism which produces this local rotation is viscosity, or fluid friction. 
Viscosity is a feature of a fluid which can be measured, and may be thought of 
as describing the ‘degree of stickiness’ of the fluid. 

Clearly, some fluids are more viscous (‘sticky’) than others, as exemplified by 
the different behaviours of treacle and water at normal room temperatures. 
However, all real fluids have some viscosity, and an important consequence of 
this is that all real fluids adhere to the boundary of any obstacle past which they 
flow. 

Thus an improved model, which takes viscosity into account, must impose a 
stronger condition on the behaviour of the flow velocity close to the boundary. 
Roughly speaking, what we now require is that the boundary should ‘consist 
entirely of stagnation points’. There are indeed more advanced models of fluid 
flow which include the effects of viscosity. However, since the locally 
circulation-free condition must be dropped, these models cannot use complex 
analysis in the way that has been done here, as pointed out in the remark 
following Problem 2.1. 

All is not lost, however, as regards the application of our first model. Often the 
effects of fluid viscosity are significant only in a narrow region called the 
boundary layer, in the immediate vicinity of the obstacle. It is then reasonable 
to ignore the viscosity of the fluid except within the boundary layer, and to 
apply the original model to flow in the exterior region. Appropriate matching 
of models is needed at the interface between these two flow regimes. For 
example, the successful prediction for an aerofoil of increasing lift for increasing 
(but small) angle of attack corresponds to cases where the boundary layer 
around the aerofoil is narrow. (In fact, the aerofoil shape is of interest and 
practical use primarily because it achieves this narrow boundary layer.) At 
larger angles of attack, the boundary layer separates from the upper surface of 
the aerofoil (see Figure 4.6), leading to a large area in which the predictions of 
an inviscid (non-viscous) model are invalid. 



The circulation which is present for an aerofoil, and creates the lift on it, arises 
from the action of viscosity within the boundary layer. Similarly, the 
circulation considered in the context of a circular cylinder (with a disc as 
cross-section), in Subsection 3.1, may be thought of as created by the cylinder 
being spun in the clockwise direction (for negative c), with the ‘obstacle 
boundary’ now being at the outside of a narrow boundary layer. 


Viscosity is caused by 
intermolecular forces acting 
within the fluid. 


This requirement, which 
amounts to taking the limit of 
flow velocity to be zero cis each 
boundary point is approached, 
is called the no-slip condition 
for viscous fluids. 
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The lift force predicted by our model in such a case is certainly observable, and 
its presence is known as the Magnus Effect. The direction of this force is 
always from the side of the obstacle where the circulation and uniform stream 
oppose each other, and towards that where they reinforce each other. 

For example, if a football is kicked from right to left, say, with a spin which is 
clockwise when viewed from above (see Figure 4.7), then its trajectory will 
deviate to the right of its initial direction, producing the so-called ‘banana 
kick’. The effects of top-spin or back-spin on tennis balls may be predicted in a 
similar manner. 


trajectory 

ball 


initial direction 



ball with 
clockwise 


Figure 4-7 


The Magnus Effect was at one time put to use in a form of ship propulsion, 
using tall vertical cylinders mounted on deck and known as Flettner Rotors. In 
the presence of a wind, rotation of the cylinders produced a horizontal ‘lift’ 
force, at right angles to the direction of the wind. Rotating cyhnders have also 
been attached to ship rudders, providing dramatic improvements in 
manoeuvrability at low speeds. 

We have looked so far at the shortcomings associated with our modelling 
assumption that the flow is locally circulation-free. We also assumed that the 
flow is locally flux-free, and pointed out in Subsection 1.2 that, for a fluid of 
constant density, this condition is equivalent to the Principle of Conservation of 
Mass. 


Heinrich Gustav Magnus 
(1802-70) was a German 
chemist and physicist, and a 
Professor at Berlin University. 
He explained the Magnus Effect 
in 1852, after interest in the 
subject had been aroused by 
deviations to the paths of 
bullets in a cross-wind. 

A football does not have 
constant two-dimensional 
cross-section, but the same 
effect is observed for spinning 
objects which do. As far as the 
ball is concerned, the oncoming 
airstream moves initially from 
left to right. The ‘lift’ for the 
football here is actually in the 
horizontal direction. Top-spin 
produces a downward ‘lift’ force 
and back-spin an upward one. 

Anton Flettner was a German 
naval engineer, and built a ship 
with these rotors in 1924. 
Experiments with the ship were 
satisfactory, but ships propelled 
by ordinary engines proved to 
be economically superior. 


A fluid of constant density is said to be incompressible. In reality, no fluid is 
completely incompressible (it is compressibility which permits sound waves to 
travel through a fluid), but liquids are much less compressible than gases. Since 
fluids are compressible to some degree, they axe not locally flux-free in a precise 
sense. 


Where aerofoils axe concerned, the predictions of models based on the 
assumption of incompressible (locally flux-free) flow axe reasonable for 
airstream velocities below 100 ms“^, which is much less than the speed of 
sound (about 330 ms“^). At higher velocities, the theory of compressible flows 
is required. The compressibility of air has the effect of increasing slightly the 
lift for a given angle of attack, and of raising the extent to which the aerofoil 
can be tilted before stalling occurs. 

Mathematical modelling is an iterative process. It is in the nature of modelling 
that a first model will be only partially effective in its task, and that 
subsequent improvements will probably be required. This implies no 
devaluation of the first model, since it is necessary to make sufficient 
simplifications when starting to consider a situation to ensmre that some 
progress is made. Within its hmitations, the model of fluid flows introduced in 
this unit may be judged relatively successful, since it predicts significant 
physical phenomena which are observable in practice and allows us to see the 
direction in which a further model must improve on the original. 





4.4 Derivation of Bernoulli’s Equation 

We assume, as earlier in the unit, that the flow is two-dimensional and steady, 
and that the fluid has constant density. As you may recall, Newton’s Second 
Law in single-particle mechanics states that 

force = mass x acceleration. 


Here we require a version of this law which is applicable at each point of the 
fluid. This can be obtained by considering both the force acting on and the 
mass of a small volume of fluid, and then letting the volume size tend to zero. 
This gives, at each point, 

force per unit volume = density x acceleration, 
or in symbols, 

(4.7) 

where F{z) is the force per unit volume at z and p is the constant fluid density. 

The velocity of the fluid is given by dz/dt = q{z), so that its acceleration may 
be expressed as 


d^z _ dq{z) 

~ ~dr' 


(4.8) 


Using Equation (4.8) to substitute for d?z/dt^ in Equation (4.7), and then 
taking the complex conjugate, we have 




^ dt 


Since q is an analytic function, this can be expressed as 




or simply as 


F{z) ^ pq{z)q'{z). 


(4.9) 


This expression for the conjugate force per unit volume acting on the fluid will 
be of use shortly. 


We now relate the force per unit volume acting within the fluid to the pressure 
of the fluid. We shall ignore the effects of ajiy other contributions to the force 
per unit volume, such as those due to gravity or to viscosity within the fluid. 

Suppose that 7(5) is a unit-speed parametrization of any line segment within 
the flow region of the complex plane, and consider a small cylinder, of 
cross-sectional area A, whose axis lies along the line segment. The ends of the 
cylinder are at 7(3) and 7(3 + h), as shown in Figure 4.8, so that the length of 
the cylinder is h. 

The force per unit volume F{‘y{s)) leads to an approximate force F{'y{s))Ah on 
the cylinder of fluid, where the relative error of the approximation decreases as 
A and h tend to zero. We shall restrict attention to the component of this force 
in the direction of 7'(s), which is Ft(s)A/i. 

The fluid pressure causes a ‘force distribution’ acting at each point of the 
cylinder in a direction perpendicular to the cylinder’s surface. However, the 
directions of the forces acting on the curved part of the surface axe 
perpendicular to the direction 7'(s), and hence do not contribute to the force 
component in this direction. The remaining forces due to the pressure are 
approximately p{‘y{s))A, in the direction of j'(s), and p(j(s + h))A, in the 
opposite direction, as shown in Figure 4.9. 


This subsection is intended for 
reading only. 


Note that F(z) is the force 
acting at the point z on the 
fluid rather than on an obstacle. 



Figure 4.8 


ph{s + h))A 



Ph{s))A 
Figure 4.9 
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It follows that 


FT{s)Ah S p('y{s))A - p('y(s + h))A; 
that is 

^^(3) ^ p(7(3))-p(7(^ + h)) 

h 

The relative error implicit in use of the ^ sign tends to zero as h tends to zero, 
and so in this limit we obtain 

=-(4-10) 

Now 7 '(s) has modulus 1 so, by the first of Equations (1.2), the component 
Ft{s) of F('y{s)) in the direction of 7 '(s) may be written as 

Fr(s)=Re(F(7(s))7'(5)). (4.11) 

On applying Equations (4.11) and (4.9) to Equation (4.10), we find that 

Re(pg(7(s))q'(7(s))7'(s)) = 
that is. 


This can be rewritten as 

= |P^(9(7(s))9(7(s))). 


Equations (1.2) were expressed 
in terms of the fluid velocity 
function q, but they apply also 
to any other function. 


Note that 

Re(ziZ2) = 5(ziZ2 +5122). 


By writing q{z)q(z) = \q{z)\^, and integrating with respect to s, we obtain 
Bernoulh’s Equation, 


P(7(s))+ ip|g(7(s))|^ =Po, 

where po is a constant. Since our derivation referred to an arbitrary line 
segment in the fluid, Bernoulli’s Equation holds throughout the fluid; that is, 

pW + |pk(2)l^ =Po- 


EXERCISES 


Section 1 

EX6rciS6 1.1 Find the component of q{z) = (at any point z) in the 

direction specified by each of the following complex numbers. 

(a) e-2‘V3 (b) _2e-W3 

EX6rciS6 1.2 Use Equations (1.2) to show that, if 6 is any angle, then a flow 
velocity q{z) can be expressed in terms of its components by 

<li^) = i^eiz) - iq{e-^/2){z))e'^■ 
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Ex@rciS6 1.3 Consider the flow velocity function 


(a) Let r be the unit circle {z : \z\ = 1}. Use Equation (1.7) to show that q has 
zero circulation along F, and zero flux across F. 

(b) Show that 5 is a model flow velocity function, using in turn: 

(i) Theorem 1.2; (ii) Theorem 1.3. 


Section 2 

Exercise 2.1 

(a) Draw a rough sketch of the streamlines for a uniform flow with complex 
potential 2z, on which is superposed a source of strength 47r at the origin. 
Find the position of the single stagnation point of the flow. 

(b) The curved streamlines which approach the stagnation point may be taken 
as the boundary of a blunt object placed in the uniform flow. Find 
equations to describe this boundary, and calculate the ‘maximum width’ of 
the object perpendicular to the flow. 

Exercise 2.2 Prove that, if a model flow velocity function with domain C is 

bounded, then the flow is uniform. 

{Hint: Use Liouville’s Theorem.) 


Section 3 

Exercise 3.1 By considering the stream function of the flow with complex 
potential 

n(z) = z + y (z&C-Ka), 

where Ka = {z : |z| < a}, find the equation of the streamhne which passes 
through the point 2ai. Hence find the limiting value of the distance between 
this streamline and the real axis, far from the disc Ka- 

Exercise 3.2 Let Ka denote the closed disc {z : \z\ < a}. Consider the 
obstacle (see the figure in the margin) 

iF = [-2,-l]Ui('iU[l,2]. 

(a) Show that the Joukowski function J is a one-one conformal mapping from 
C-K onto C-[-§,§]. 

(b) Find a value of a for which the function 

Ja(u;) = w + — 

is a one-one conformal mapping of C — Ka onto C — [— |, |]. 

(c) Show that the Obstacle Problem for K, with circulation 27rc around K, has 
solution 

“ J{zWl-25/i2J{zw) 

{Hint: You may assume that the function / = o J satisfies 
Condition (3.5). Also, find the velocity function from the complex 
potential function, using the fact that Da.c = Ja — ic Log and the approach 
of Problem 3.5.) 


Unit B2, Theorem 2.2 
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ExerciSG 3.3 Suppose that 0 < a < r. Let K be the obstacle whose boundary 
is the ellipse 


r = {. 


(r-ayrf (r + «Vr)' 
whose major axis Ues along the imaginary axis (see the figure in the margin). 


.-)■ 


Show that the Obstacle Problem for K, with circulation 27rc around K, has 
solution 


V + J 


{Hint. Find the flow velocity function from the complex potential function.) 



Section 4 

Exercise 4.1 in Problem 3.7, you found that the model flow past the ‘plate’ 
K = [—2a, 2o], with circulation 27rc around K, has a flow velocity 

Find the force on the plate by applying Equation (4.4), the modified form of 
Blasius’ Theorem. 

{Hint. The circulation around K is given by 



where F is any simple-closed contour surrounding K.) 

Exercise 4.2 The asymmetric Joukowski aerofoil K, shown on the left of the 
figure below, has boundary dK = J{dD), where 

The condition Re < 0 ensures 
that the critical point w = —1 
lies within D, with its image 
z = -2 within K. For Re/? = 0 
(but Im/3 > 0), the circle dD 
would pass through both of the 
critical points w = ±\. In this 
C2ise, the image J{dD) is a 
circular arc in the upper 
r-plane, from r = —2 to r = 2, 
and J is a one-one mapping of 
— D onto C — J{dD). The 
case for Re/3 = 0 is indicated by 
the broken curves on the 
diagram. 


D = {w.\w — P\ <\1 — PW and Re/?<0. 

It can be shown that the Joukowski function J is a one-one mapping of C — D 
onto C — Ff. 



Show that, if the Kutta-Joukowski Hypothesis holds, then the force acting on 
K due to the flow is 

F = dTrpi Im /8. 
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SOLUTIONS TO THE PROBLEMS 


Section 1 


1.1 (a) For q{z) = z, all of the arrows for z ^ 0 have 
direction Arg z, and so point radially outwards. At the 
origin, z = 0, there is a stagnation point. 




(b) By using results from Unit Al, Subsection 2.3, we 
have that the direction of q{z) = i/z {z ^0) is given by 
Arg(i/z) = Argi - Argz + 2Ti7r 
= Argi + Argz + 2n7r 
= Arg z + 7r/2 + 2Tiir, 

where n € Z. The pattern is therefore as shown below. 
This is an example of a vortex flow. 


J 

\ N 
t t t . 

nr 

-* / f 

\ 





1.2 In each case, the streamline picture can be deduced 
from the corresponding arrow pattern in the solution to 
Problem 1.1, as shown in the following figures. 


(a) 



(b) 



1.3 With q{z) = 5e and 6 = 27r/3, we have, from 
Equation (1.1), 

q2ir/3{z) = 

= Re(5e®’'/®) 



= -|V3. 


1.4 Prom Equation (1.1), we have 
qe{z) = Re{^e*0). 

With 6 — 1^12 in place of 6, this becomes 
q(e-^/2){z) = Re(g(z)e‘^*“”^*>) 

= Re(^e‘®e""/^) 

= R«(-ig(z)e‘®). 

Now if w = u + iv is any complex number, then 
Re(-iu;) = Re(v - iu) = v = Imiy. 

It follows that 

9(9-V2)(«) = Im(g(z)e‘®). 


1.5 (a) If q{z) = z is the conjugate velocity function, 
then the velocity function itself is ^(z) = z. To draw the 
arrow pattern for this velocity function, note that z has 
the same real part as z, but an imaginary part of 
opposite sign. This gives a sketch like that below, using 
arrows of fixed length (except at 0, which is a stagnation 
point of the fiow). 




/ 

\ 

:: 


4 = , , ' 


4 = ! 


t 

- 

\ 

/ 




(b) If q{z) = 1/z, then we have 
9(z) = i/z (^ec-{0}). 

Arg(9(z)) = Arg(l/z) = Argz, 
so the direction of g(z) is radially outwards. This leads to 
the same fixed length sirrow pattern as in the solution to 
Problem 1.1(a), except that the origin should in this case 
be omitted. 

(c) If g(z) = i/z, then we have 
q(z) = -ilz (z€C-{0}). 

The picture is like that of the solution to Problem 1.1(b), 
but with all the arrows reversed. 
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1.6 (a) The velocity function q{z) = ^ (z € C) has 
conjugate function q{z) = z. This is not an analytic 
function (see Problem 1.8 of Unit A4), so, by 
Theorem 1.2, q is not a model flow velocity function. 

(b) For q{z) = ijz (z € C — {0}), the conjugate velocity 
function is q{z) = —ijz. This function is anailytic on 
C - {0} (see Problem 1.2 of Unit A4, and apply the 
Multiple Rule in Theorem 1.3 of that unit), so, by 
Theorem 1.2, g is a model flow velocity function. 


1.7 (a) The flux of q(z) = Xjz across the unit circle T is 
(from Equation (1.7)) 

.Fr = Im J q{z) dz 
= Im y* “ 

= Im(27ri) (by Cauchy’s Integral Formula) 

= 27r. 

(b) Similarly, the circulation along F of q{z) = —f/z is 
Cr = Re J q{z) dz 

= Re J - dz 

= Re(—2x) 

= -27r. 

(c) In each of parts (a) and (b), the origin is excluded 
from the domain TZ of q, so that the inside of F (the unit 
circle) does not lie in Ti. The function q is, therefore, still 
locally flux-free and circulation-free, in each case. 


1.8 (a) By the Cauchy-Riemann Theorems ( Unit A4, 
Subsection 2.1), if the partial derivatives of the real 
functions u, v with respect to x, y exist and are 
continuous on TZ, then u -f- in is an analytic function on TZ 
if and only if 


du dv 


du 


n7^. 


Now the velocity function g = gi -t- ig2 has conjugate 
function q = qi — iqz, and g is a model flow velocity 
function on TZ if and only if g is an2dytic on TZ (by 
Theorem 1.2). Putting these results together, with u = qi 
and V = —g2, we have that g is a model flow velocity 
function on TZ if and only if 

^ + ^=0 and ^-^=0 onTZ. 

dx dy dx dy 

(b) The real and imaginary parts of 

. ^ + iy) 

■ |zP x^ + y^ ’ 
i, respectively. 


g(z) = 


gi{x, y) = - 


92 (x,y) = - 


The partial derivatives of gi and g2 are 

dqi _ 2xy dqi _ y^ - x^ 

dx (a;2-f.j,2)2’ 9y “ (x2 + j,2)2 ’ 

dq 2 _ y^ — x^ dq 2 _ 2xy 

dx~{x^+y 2 f’ dy ~ (x2+y2)2- 

Hence we And, if z ^ 0, that 

$1 + $1 = 0 „d = 

ox dy dx dy 


By the result of part (a), g(z) = i/z (z e C - {0}) is 
therefore a model flow velocity function. 


1.9 (a) For the velocity function g(z) = z, we have 

9i(X) y) = X and g2(x, y) = y, so that 
dqz _ _ 

dx dy 

Hence, by Theorem 1.3, g is locally circulation-free. 

However, since 

the function g is not locally flux-free. 

(b) For the velocity function g(z) = iz, we have 

gi{x,y) = -y and q 2 (x,y) = x, so that 
dqi dq 2 
dx ^ dy 

showing g to be locally flux-free, by Theorem 1.3. 

However, since 
dq2 _ dqi _ 
dy ' 

the function g is not locally circulation-free. 




Section 2 

2.1 The complex potential Q is deflned by 
n'(z) = g(z). 

With n(z) = ^{x,y) + i^{x,y), we have 


dx 


(by the Cauchy-Riemann Converse Theorem) 
~ dy ^ dx 

(by the Cauchy-Riemann Theorem). 

On taking complex conjugates, we obtain 


g(z) 


a’F .d<ff 


dy 


dx ' 


2.2 A complex potential Q is given by n'(z) = g(z). 
Here we have 

9(^)=zi^ (z€C-{ 0}), 
so we require for Cl a primitive of 



We need to choose a simply-connected subregion of 
C — {0} to be the domain of fl, as discussed on the audio 
tape (see Frame 3). If the domain is chosen to be the cut 
plane C,, we may take 
n(z) = -(1 -|-8i)Logz 

= -(1 + 8i)(log. |z| -I- i Arg z). 

The stream function is then 

Im(n(z)) =-(81og. |z|+Argz). 

The equation for streamlines is Im(f2(z)) = k, where k is 
a constant; that is, 

-81og^ |z| = A: -I- Argz. 

Further manipulation of this equation gives 
|z| = 

where K = is another constant. As Arg z increases 

from —TT to TT, |z| decre2ises exponentially. On joining this 
picture up across the negative real cixis (using a 
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logarithm function defined on a different cut plane), we 
see that the streamlines are spirals directed inwards, as 
shown below. This is a whirlpool or spiral vortex flow 
(combination of a sink and a vortex). 



The source strength and vortex strength are, respectively, 
the flux across and circulation along any simple-closed 
contour F which surrounds the origin. These Me, in turn, 
the respective imaginMy and real parts of 

j q{z)dz = j 

which, by Cauchy’s Integral Formula, has the value 
—2iri(l -h 8i) = IAtt — 2i:i. 

Hence the flow has source strength —27r and vortex 
strength IAtt. 

2.3 (a) We have to show that 
= q(z), for 2 € 7^, 
where 72. = C — {a: € K : x < /i}. On 72, 
n'{z) = Log'( 2 ) - Log '(2 - h) 

_ 

~z z-h’ 

so that 

as required. 

(b) The complex potential 
Q{z) = Log 2 - Log (2 - h) 
leads to the stream function 

Im(n( 2 )) = Arg 2 - Arg (2 - h), 
which must be constant along each streamline. Suppose 
first that 2 is in the upper half-plane. 



The angle <!> marked on the diagram above is equal to 
Arg(2 — h) — Arg2, since the exterior angle of any 
triangle is equal to the sum of its two interior opposite 
angles. Hence as z moves along a streamline, the angle 
remains constant. It follows that z moves along a circulM 
arc paissing through 0 and h. By symmetry, the centre of 
the circle lies on the line x = \h. 


A similM conclusion is reached by considering z to be 
below the real aixis. The real axis itself contains three 
streamlines, namely ]-oo,0[, ]0,/i[ and ]/i, oo[. The 
directions on the streamlines follow from the facts that 
2 = 0 is a source and 2 = /i is a sink. The streamline 
diagram is therefore as follows. 



2.4 Taking the limit, we have 

_ _J_ 

which is the velocity function for the doublet. The 
streamlines for the doublet Me shown in Frame A of the 
audio tape. 


Section 3 

3.1 (a) The condition for stagnation points is q{z) = 0, 
or equivalently, q{z) = 0. Here we have 

9(2) = !-^, for2€C-{0}, 

so that the stagnation points satisfy 
z^-a^ = {z + a){z-a) = 0. 

Hence the only stagnation points are 2 = ±0. 

(b) Substitution of 2 = ae‘‘ into q{z) = 1 — 0^/2^ gives 



= (-2isint)e“, for 0 < t < 27r. 

This gives the velocity of the flow at 2 = ae’‘, on the 
circle {2 : |2| = a}. The flow speed here is 
|g(ae“)| =2|sinf|. 




3.2 The quadratic equation 


2=5 ^ic ± y/Aa^ — , 

ich are the stagnation points of the flow. 


(a) For — 2o < c < 0, we have 4a^ — > 0, so that the 

stagnation points have real parts ±5^40^ — and 
imaginaxy part 5C. The modulus in each case is 
|2| = iV'(4a=*-c2) + c2=a, 
so that both stagnation points lie on the circle 
{z : jzl = a}. They are below the real axis (since c < 0), 
and symmetric2Jly placed on either side of the imaginaxy 


(b) For c = —2a, we have 4o^ — = 0, and so 

z = jic = —ia. Thus there is a single stagnation point, 
which lies at the intersection of the circle {z : \z\ = 0} and 
the imaginary axis. 

(c) For c < -2a, we have 4a“ - c^ < 0, so that the 
stagnation points axe at 

z = \i{c ± y/— 4a*) = 2±, say. 

Both of these stagnation points are on the imaginary 
axis, and below the real axis (since c < 0 and 
|c| > y/c^ — 4a^). Now since c < —2a, we have 
Imz- = Kc - y/(? - 4a2) < 5C < -a, 
so that [z-l > a, and z_ lies outside the circle 
{2 : \z\ = a}. Also 2+2_ = —a^, so that 
12+1 = a^/|z_| < a, and z+ lies inside the circle 
{2:|2l=a}. 


3.3 We show that q{z) — 1 satisfies the three properties 
given in the statement of the Obstacle Problem. 

Property (a) Since q{z) = 1, lim q{z) = 1. 

Property (b) A complex potenti2J function for q is 

n(2) = 2 (zeC-K). 

Let aedK = K. Then 
hm Im(n(2)) = k, 

where k = lm a, since Im Q h2is continuous extension to 
the whole of C. 

Note that Im a is constant on dK because dK is a line 
segment parallel to the real axis. 

Property (c) Since q has analytic extension to the whole 
of C, we have Cr = 0, for any simple-closed contour F 
surrounding K. 


3.4 The translation /(z) = 2 — /3 is a one-one conformal 
mapping from C - K onto C- Ka- Also / satisfies 
Condition (3.5), with ao = —0 and a_i = a_2 = ... = 0. 
Thus, by the Flow Mapping Theorem, the velocity 
function which solves the Obstacle Problem for K, with 
circulation 27rc axound K, is 

= qaAf{z))TF) 


'{Z-PY Z-p' 


3.5 By the Flow Mapping Theorem, the required 
complex potential function is 
n(2) = na.c(/(2)), 

where /(z) = J-^{z) = §(2 -t- z-^l - 4/2*) and a = 1. 
Using the hint, we have 
n(2)=ni.e(j-‘(2)) 

= J(J“‘(2)) - icLog(j"‘(z)) 

= 2 - icLog(i(2 -I- 2^1 - 4/22)) , 


as required. 

The velocity function is q{z) = n'(2), where 
ic(j-i)'(2) 

J-i(2) 

ic(l - l/(J-‘(2))' 


9 !{z) = 


J-Hz) 


(as in Example 3.1) 


J-i( 2)-1 /J-Hz)- 

Hence, as in Example 3.1, 




=■-( “ y 

\zyfrr^) 


3.6 (a) The sc2iling 21 = zja maps C — Ka onto 
— Ki \ the Joukowski function 22 = J{i\) = ^1-1- l/«i 
maps — K\ onto C — [—2,2]; the scaling w = azi maps 
C - [-2,2] onto C - [-2a, 2a]. 

Since each of the above mappings is one-one, onto and 
conform2J, the function 




^ \ 

is a one-one conformal mapping from C - o onto 
C- [-2a,2a]. 

The set dKa. is |ae**: 0 < f < 2'7r}. The image under Jo 
of a typical point in this set is 
Ja(ae“) = ae*‘ -I- a^/(ae“) 

= a(e“ -|-e~“) = 2acost. 

Hence 


Ja{dKa) = {2acost: 0 < f < 27r} 

= [-2a, 2a], 

as required. 

(b) The inverse function z = Ja^{w) is the composition 
of 

the scaling 22 = w/a; 

the inverse Joukowski function 

zi = J~^{z2) = 5(22-1-22\/1-4/2|); 
the scaling 2 = azi. 

Hence the required inverse function is 

Ja^{w) = ^(w + wy/l^^ 4 c^Jw^) (w 6 C — [—2a,2a]). 
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3.7 We know, from Problem 3.6, that /(z) = Ja^{z) is a 
one-one conformal mapping from C — [—2o, 2a] onto 
C — Ka- Now we show that / = J~^ satisfies 
Condition (3.5). 


The binomial series for (l - is 

(i-VA-)'" 

for |—4a^/z^| < 1, 

= 1 — ^^— • • •, for |z| > 2a. 





= Z-- -^-, for jzj > 2a, 

so that Condition (3.5) is satisfied. 


It follows from the Flow Mapping Theorem that the flow 
velocity function is 


q{z) = qa4J:\z)){Ja^y{z) 


By the Inverse Function Rule, 

*(i-“V(J.-W)’)" 

{since J!,{w) = l- a^/w^) 

SO that 


^ {jz\z)-ayjr\z)) 

=i-( . " ) 

\z^l-4ayz^J 


(since ^(z) + a^/Ja ^(^) = z). 


3.8 (a) The given circle is 

Cr = {w: |ii;| = r} = {re“ : 0 < f < 27r} . 
The image under Ja of a typical point in Cr is 



cost-b 




The image set Ja(Cr) is therefore 




sint: 0 < f < 27r 



This is an elhpse in the z-plane, with semi-major axis 
r -b a*/r and semi-minor axis r — a^/r, as shown in 
Figure 3.9. 

(b) The function is a one-one conformal mapping 
from C - [-2a, 2a] onto K^, where Ka = {w: |u;| < a}. It 
maps the boundary dK = JaiCr) of the given obstacle to 
the boundary dKr = Ct of the disc Kt = {w: |u;| < r}. 
Hence the restriction / of to C - R is a one-one 
conformal mapping from C - K onto C- Kt- 
Also, by Problem 3.7, we know that / satisfies 
Condition (3.5). 

It follows from the Flow Mapping Theorem, and using 
the Inverse FNmction Rule as in the solution to 
Problem 3.7, that the flow velocity function is 

q{z) = qr4JZ\z)){Ja^y{z) 


/ {Jr\z)y -r^- icJa^z) \ ( {Jl\z)y \ 

V )\{J-a\z)y-a^) 


( {Ja. ^{z)y -a^ +a^ -r^ - icJa '(z) 

1 , {Ja-Hz)y-a^ , 


^ / r^-a^-bicJa^(z) \ 

' 1 , {Ja-Hz)y-a^ )’ 

where J^^iz) = \{z -b z-y/l - Aa^jz^). 


3.9 (a) The function Ja ^ maps the boundary dK of 
the aerofoil onto C, as well as mapping - K one-one 
onto C — D, where D is the closed disc with boundary C. 
Furthermore, the translation g{w) = w + b maps C to 
dK(a-\-b), where R(a+6) = {in : lto| < a -b 6}, as well as 
mapping £, — D one-one onto C — iC(a+6)- Thus the 
composite function 

nz)={goJ:^){z) = JZ\z) + h 
maps dK to dK(a+b) , and also maps — K one-one onto 
C-K^a+b). 

Since both and g satisfy Condition (3.5), so does /. 
By the Flow Mapping Theorem (with a -b ft in place of a), 
a suitable complex potential for the flow described is 
n(z) = (n«.+i),,o/)(z) 

= n(a+b)4Ja\z)-^h) 

= w -b 6 -b -ic Log(w -b 6), 

w = Ja^{z) = \{z -b z-^/l - Aa^/z^). 

(b) Using the Chain Rule (or directly from the formula 
given in the theorem), we obtain 

q{z) = (^1 - 

where, as in Problem 3.7, 
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It follows that, in terms oi w = ^(z), the flow velocity 

function is 


9(^) = - 


(. _ (a + b)^ 

{w + by u 
e have w = Ja^{z) 


.)■ 


9(^) = - 


/ (w + g + 2b){w - g) 

I, {w + by 


f w^w + a + 2b) icw^ \ 

\iw + a){w + by {w + b){w^-a^)J- 
Thus if c 7^ 0, then lim q(z) does not exist. If c = 0, 
then the limiting velocity aa 2 —» 2g is g/(g + 6). 


3.10 (a) Prom Equation (3.10), since a = ae and 
g > 0, we have 

Ja = R^^oJ^oR^. 

If 5 is the circle R^^{C), then 
C = R^{B). 

It follows that 

dK = {R^^ o J„) (C) 

= {R^^oJ^oR^){B) 


where 

P = 2{a + - ic, 

Q = (g + 6)((g + - l) - ice-**) . 

Now the limit 

lim q{z) = lim q{z) = lim^g(2) 
exists if and only if Q = 0, that is, if and only if 
(g + 6)(e-=‘»‘^-l) 

= -t(g + 6) - e'*) 

= -2(g + 6)sin^. 

(For this value of c, we have 


lim q{z) = 


a^P 

2a{a + bye-^'* 


ae** (2(g + b)e-** + 2i{a + 6) sin (f>) 
2 (g + by 


The complex conjugate of this is 

, , ae-'*cosS a cos 6 

a+b =inT- 

As <t>—*0, this expression tends to the value g/(g + b) 
found in the solution to Problem 3.9.) 


(b) Let D be the closed disc with boundary B which has 
centre —be-'* and raxlius g + 6. The function ^ maps 
C — K onto C — D, and the translation g{w) = w + be-'* 
maps — D onto C — K^a+b)- 
The restriction / of the function 
(5oJ-i)(2) = J-‘(2) + 6e-‘* 
to C — if is a one-one conformal mapping from C — if 
onto C — if(0+6). 

Also, / satisfies Condition (3.5) since g and do. 
Hence, by the Flow Mapping Theorem, a suitable 
complex potential for flow past if, with circulation 27rc 
around if, is given by 

n(2) = (n(„+^).,o/)(2) 

= ii(a+b),c{J^\z) + be-**) 

= w + be-'* -b - tcLog(w + be-'*) , 

where w = J-^{z). 


(c) The method of Problem 3.7 can be used to show that 

Following the hint, we use the variables 


W2 = w + be '* = wi + (a + b)e '*. 

We may then write, by the Flow Mapping Theorem, 




u)i(tt;-|-Q!) 


w* (w| - (g + by - icw2) 
wi{w -f a)w2 

(wl + Pwi + Q) 

Wi {w + q)W2 ’ 


Section 4 

4.1 If the ball is pushed sideways from the centre of the 
airstream, then the flow speed on the side of the bedl 
closest to the stre2un centre will be greater than that on 
the opposite side, since the latter part of the flow 
encounters more resistance to its motion due to the ball. 
Bernoulli’s Equation predicts that there is, in 
consequence, a higher pressure on the side further from 
the stre2im centre, amd the pressure difference provides a 
force (from the higher pressure region to the lower) which 
tends to restore the ball to its central position in the 
airstream. 

4.2 (a) By applying Bernoulli’s Equation, 

P(«) + 5Pk(«)l* =Po, 

to the velocity function 



we obtain 

p(^)=P0-f(l-f2-7) + 

(b) On putting z = ae'* and z = ge~’* into the 
expression for p{z) in pcirt (a), we find that 

Piae**) = Po - ^ “ ^^^) 

. p. - e 2) (.-«-.« + ^) 

= Po -b |(2tsint- 
= Po - |(2sint-. 
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(c) The symmetry of p{ae'*') about the imaginary axis 
follows from the fact that 
sint = sin(7r-t). 

If c = 0, then there is symmetry of p(ae’‘) also about the 
real axis, because 
sint = -sin(-t), 
so that 

sin*t = sin*(—t). 

(These symmetry results also apply to the flow as a 
whole. We have 

viz) = Pit), 

for symmetry about the real axis, when c = 0, and 
viz) = Vi-z), 

for symmetry about the imaginary axis.) 

4.3 The required total force is 

F = —a J ^po — ^ ^2sinf ~ ^ e'*dt. 

Now we have 

-“I 

so the real part of F is 

ReF = 2pJ (asin^t - csint) costdt 

= 2p[iasin®t — icsin^t] 

= 0 , 

and the imaginary part of F is 

ImF = 2p J (asin® t - csin® f) dt 

= J (“(l “ *) "• 2*^(1 ~ cos 2t)) dt 

— 2p[—a(cost — 5 cos® t) — ^c{t — i sin2t)]^’^ 

= —27rpc. 

Hence we have F = —2‘Kpci. The lift is —2npc (positive if 
c < 0), and the drag is zero. 


4.4 Applying Blasius’ Theorem, the conjugate force is 





= 27rpci (by the Residue Theorem). 

The total force on the cylinder is therefore F = -27rpci 
(per unit length into the page). 

4.5 The lift is given by the Kutta-Joukowski Lift 
Theorem as —27rcp. Also, from Problem 3.10(c), the 
Kutta-Joukowski Hypothesis holds if and only if 
c = —2(a + 6)sin<^. 

In this case, the lift on the aerofoil is 

-2irp{-2{a + b) sin<^) = 47rp(a + 6) sin <i>. 

This lift is positive for 0 < (^ < jr, and increases with 4) 
for 0 < <^ < ir/2. 

In practice, this prediction for the lift is borne out only 
for small values of the angle of attack tj). Above 15° or so, 
there is a rapid decrease in the lift provided (stalling), 
and flight can no longer be sustained. 
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SOLUTIONS TO THE EXERCISES 


Section 1 


1.1 (a) From Equation (1.1), with g(z) = and 

6 = —ItcIZ, we have 

9_,./3(z) = Re(3e--/-e--/^) 


= Re(3e 



= -3/V2. 


(b) The second direction is specified by 
SO that the required component is 

= Re(3e*"/‘‘) 

= 3 COS - 
4 

= 3/^. 


1.2 We have 

q{z)e'^ ^ Re{q{z)e'^) + ilm{q{z)e*^) 

= qe{z) + iq{B-r/2){z) (by Equations (1.2)). 
Taking the complex conjugate of both sides, we obtain 
q(z)e~'^ = qe{z) - iq(e-^/2){z), 
so that 

q{z) = {qe{z) - iq(e-ir/2){z))e'^. 


1.3 (a) The conjugate velocity function is 

Qi^) = 1 - 

so that, from Equation (1.7), we obtain 
Cr + i:Fr = j (l - dz. 

This integral can be shown to be zero either by using the 
standard parametrization for the contour T, or by 
applying the Residue Theorem. Hence we have 
Cr = .Fr = 0. 


(b) (i) The conjugate velocity function q is analytic on 
C — {0}. Hence, by Theorem 1.2, g is a model flow 
velocity function on this region. 


(ii) The real and imaginary parts of 

\* (x2+j/2)2 


are, respectively, 
qi{x,y) = 1 - 


92(x,j/) = 


2xy 


{x^+y^r (a:2 + y2)^- 

The partial derivatives of qi and 92 are 

dqT, _ 2x{x^-3y^) dqi _ 2y{3x^ - y^) 

9x (a;2+j^2)* ’ dy (a:2+j/2)3 

dq2 _ 2y(3x^ - y^) dq2 _ 2x(3y^ - a;^) 
dx (a;2^j^2)3 ’ dy (a;2^.y2)3 

Hence we find, if z 5^ 0, that 

dx'^ ^ 


By Theorem 1.3, q is locally both flux-free and 
circulation-free; that is, 9 is a model flow velocity 
function. 


Section 2 


2.1 (a) The sketch below was obtained by considering 
the flow expected when a source is placed in a uniform 
flow. (It may be confirmed by using the stream function 
derived in part (b).) 



The source of strength 47r at the origin has flow velocity 
function 2/2, so that the overall velocity function here is 

9(^) = 2(l+i). 

The stagnation point is the solution of the equation 
q{z) = 0, which is 2 = —1. 

(b) A complex potential for this flow is 
n(2) = 2(2 + Log2) (2€C0, 
and so the corresponding stre2im function is 
Imn(2) = 2(Im2-t-Arg2). 

The stagnation point at 2 = — 1 lies outside the domain 
of n, but we have 

lim^(Im2) = 0 and lim (Arg 2) = ±7r, 

where the plus or minus sign is chosen depending on 
whether the limit is taken from above or below the real 
eixis, respectively. Hence the curved streamlines which 
emerge from the stagnation point have equations 
y + Arg2 = ±7r. 

The width of the object outlined by these streamlines, 
corresponding to the point z = x + iy on the upper 
streamline, is 

2j/ = 2(7r-Argz). 

The maximum width (never att2uned, but equal to the 
limiting width as 2 —> 00 along the streamline and hence 
Arg 2 —► 0) is therefore 27r. 
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2.2 If g is a model flow velocity function which has 
domain C and is bounded, then there exists a real 
number K such that 
|9(^)|<i<:, forzeC. 

But q is an analytic function with the same domain as q, 
and is therefore entire. Also, we have 
|9(z)l = |9(z)|, for^eC, 
so that q is also bounded. It follows from Liouville’s 
Theorem that 5 is a constant function, and hence that q 
is constant. The corresponding flow is therefore uniform. 


Section 3 

3.1 The stream function is 

which is constant on each streamline. Thus the 
streamline through 2ai (at which point, x = 0 and 
y = 2a) has the equation 



As 2 -+ oo, we have (x^ + j/^) —► 0. Hence y —* §a, 

and so the distance between the streamhne through 2ai 
and the real axis tends towards the limiting value |a. 

3.2 (a) The Joukowski function J is a one-one 
conformal mapping from C,- Kx onto C — [—2,2]. Also, 
we have 

J([l,2]) = |z -t- i ; 2 e K, 1 < 2 < 2| = [2, f] , 
and, similarly, 

J([-2,-!]) = [-§,-2]. 

The required result follows. 

(b) The image under Ja of C - ifo is C - [-2a, 2a], so 

(c) The function / = o J is a one-one conformal 
mapping from C — K onto C - 1^5/4, which satisfies 
Condition (3.5), by the hint. 

Hence, by the Flow Mapping Theorem, a complex 
potential function for this Obstacle Problem is 
n(2) = (ns/4,c°/)(^) 

= J(2)-icLog(j-/,(J(2))), 
by the hint. 

Proceeding as in the solution to Problem 3.5, with 
a = 5/4 rather than 1 and J{z) in place of 2, we find that 
the velocity function is 

q{z) = I 1 - , - : I J'{z). 

3.3 The boundary ellipse here is of the same form 2is 
that considered in Problem 3.8, but rotated through an 
angle 7r/2. Hence the function 

■^-t/2 R-ir/2 = Ja, 

where a = — ai, is a one-one conform2J mapping 

from C- Kr onto C- K, where Kr = {w: |ii;] < r}. 


It follows that the function / = is a one-one 
conformal mapping from C — K onto C — Kr. Also, as in 
the solution to Problem 3.7, f = J~} satisfies 
Condition (3.5). 


Hence, by the Flow Mapping Theorem, a complex 
potential function for this Obstax:le Problem is 


n(^)=(nr.cOj-')(2) 




J:i^(z) = \(z + zy/l + 4ayz^). 

Proceeding as in the solution to Problem 3.8, we find 
that the flow velocity function is 


/ -ba^ -bfcJ^;^(2) '\ 


Section 4 


4.1 On taking the conjugate of the expression for q{z), 
and squaring, we have 




Zy/l-^a^lz^ - 4a2 


= 29(2) - 1 - - 


By applying Equation (4.4) and using the hint, we obtain 


R =\ip dz. 



where P is any simple-closed contour surrounding K. 


The remaining integral can be shown to be zero, by the 
Residue Theorem (the integrand has simple poles at ±2a, 
and the residues at these poles cancel out). Hence 
F = 27rpci, so that the force on the plate is 
F = -2npci. 


4.2 The function J~^ maps C — K onto C — D, and the 
translation g{w) = w - P maps C — D onto C — Ffii-pi, 
where R|i_p| is the closed disc with centre 0 and radius 
|1 — P\. Hence the function f = g o J~^ is a one-one 
conformal mapping from C — K onto C — K\x-g\. 

Also, since g and satisfy Condition (3.5), so does 
f = 9oJ-\ 

Hence, by the Flow Mapping Theorem, a complex 
potential function for this Obstacle Problem is 
n(2) = (n|i_^i.,o/)(z) 

= - 0 ) 

= w-P+ - icLog(tu - P), 

where w = J~^{z). 
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q{z) = 


The corresponding conjugate velocity function is 

\l-0f ic \ 

{w-PY w-p)w^-l' 
since (•/“*) («) = 'uP/{yP — l); see the solution to 
Problem 3.7, with a = 1. 

as z —► 2 (corresponding to w —> 1) 


This tends 
unless 


lim ( 1 


' {w-PY 


= 0; 


that is, unless 



= -i(l -J-{\-p)) 

= i(P-P) 


This is the value which c must take, according to the 
Kutta-Joukowski Hypothesis. By the Kutta-Joukowski 
Lift Theorem, the corresponding force on if is then 
F = —2iTcpi = Aivpi Im p. 
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